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ABSTRACT 


Acquisition, Tracking and Pointing (ATP) of space and airborne optical beams is 
fast becoming an important research topic as the requirements for the pointing and 
control of the beam is increasing. Arc-second accuracy, nano-radian jitter and large 
flexible structures require stringent pointing requirements testing the limits of control 
systems. The equipment on board the “bus” results in narrowband as well as random 
structural interactions that further complicate the control method. The effect of the 
atmosphere on the beam adds a broadband disturbance, resulting in an optical beam that 
has been corrupted by “colored noise”. A unique Laser Jitter Control Testbed was 
designed and built in order to study means to correct disturbances in a laser relay system. 
Several methods are implemented on the testbed, from classical control to adaptive 
systems, in order to develop and improve optical beam jitter control. In Adaptive control 
systems, the Least Mean Squares (LMS) and Gradient Adaptive Lattice (GAL) method 
have been used for disturbance rejection in electronic signals for years. However, these 
electronic signals are zero-mean and thus the effect of bias or phase difference in the 
error has been rarely studied. Without a compensator, a slowly varying bias in the error 
will result in a slow convergence to the steady state values. A phase difference between 
the reference signal and the error will prevent the GAL filter from removing energy in 
that tonal. A new technique called the Adaptive Bias Filter (ABF) applies a 
compensating bias to the reference signal for the adaptive filter, allowing rapid 
convergence in correcting narrowband disturbances. This technique works for both the 
LMS and the GAL filter, which were modified for use on the testbed. A second 
technique, called the Adaptive Delay Filter (ADF), was applied in a new way to the 
adaptive control systems to better correlate the error signal with the reference signal, 
resulting in faster convergence for low adaptation rates. Experiments using these 
techniques effectively demonstrated the usefulness of the new control methods, including 
the use of a combination Linear Quadratic Regulator (LQR) and LMS/ABF filter to 
remove disturbances caused by colored noise injected into an optical beam. The 
combination controller results in a 38 dB decrease in jitter caused by a 50 Hz vibration 
and a 10 dB decrease in an 87 Hz vibration of the supporting platform, as well as an 
overall 75% improvement in the broadband jitter experienced by the optical beam. 
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I. INTRODUCTION 


A. MOTIVATION 

Many future space missions, laser communication systems, and imaging systems 
will require optical beam pointing and jitter control in the sub-micro-radian regime. Jitter 
is the undesired fluctuations in the pointing of an optical beam due to the environment or 
structural interactions, and consists of both broadband and narrowband disturbances. The 
control of jitter is also a challenging problem for current programs such as the Airborne 
Laser and the James Webb Space Telescope. Jitter has a great affect on the resolution of 
an image or the intensity of an optical beam. For example, a 100 mm diameter laser 
beam with 10 //rad of jitter will result in roughly a 400 fold decrease in the intensity of 
the beam at 100 km, due to the jitter alone. In order to study this problem, the Naval Post 
Graduate School (NPS) developed the Optical Relay Mirror Lab - Laser Jitter Control 
(LJC) Testbed, located in the Spacecraft Research and Design Center (SRDC). Several 
methods are implemented on this testbed, from classical control to adaptive systems, in 
order to develop and improve optical beam jitter control. Proven adaptive means are 
tested and compared with the classical control model on disturbances provided by a 
randomly moving laser beam and vibration of the platform containing the control system. 
The control system consists of Position Sensing Detectors (PSD) and a Fast Steering 
Mirror (FSM) which are used to correct the beam’s deviations from the center of a target 
PSD. Several techniques are developed to improve the perfonnance of the adaptive 
systems. 

At the NPS SRDC there is a major research program on optical beam pointing and 
jitter control. The current application is the Bifocal Relay Mirror spacecraft. This 
spacecraft is composed of two optically coupled telescopes used to redirect the laser light 
from ground-based, aircraft-based or spacecraft-based lasers to distant points on the earth 
or in space for a variety of non-weapon force enhancement missions. Based on an Air 
Force Research Laboratory (AFRL) study, the potential missions are camouflage 
detection and penetration, chemical warfare agent detection and identification, 
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illuminator for night time/active imaging, laser fence for aircraft detection, and tunnel 
and underground structure detection. Students at NPS developed a design for the Bifocal 
Relay Mirror Spacecraft in the summer of 2000 with AFRL sponsorship. The effort 
identified several unique technologies to be developed on acquisition, tracking and 
pointing as well as beam control optics. The SRDC submitted a proposal to the National 
Reconnaissance Office to develop these unique technologies for a Bifocal Relay Mirror 
Spacecraft (BRMS) with AFRL as the co-investigator for beam control optics. The 
NPS/AFRL proposal was selected for the award in December 2000. The SRDC designed 
and had built a BRMS simulator to validate experimentally the control methods necessary 
to accurately stabilize the platform and meet the stringent optical beam pointing 
requirements. 

The SRDC has conducted research in the past on active control of narrowband 
disturbances in structures. However, the optical beam in a relay system will have 
broadband noise in addition to the narrowband jitter that results from rotating equipment 
on the spacecraft bus or aircraft platform. The object of this research is to develop 
improved control techniques to remove both the narrowband and broadband jitter in the 
optical beam. These techniques will be validated both analytically and experimentally. 

B. LITERATURE REVIEW 

1. Background 

The control of disturbances from unwanted vibration has its beginning in the area 
of acoustics. The first methods used to control vibrations were in an attempt to reduce 
unwanted noise (acoustic vibration). These methods began with the use of passive 
systems, basically blanketing the area to be protected with material that would absorb the 
frequencies of interest. Using a damping material to reduce structural vibrations in a 
similar manner is a useful technique for structural disturbances. These types of systems 
work well if size and weight are not an issue, but this is not the case for airborne or space 
borne systems. 

The earliest attempts to actively control sound are traced to Paul Lueg of 
Germany who tried to cancel sound propagating through a duct using a 180 degree phase- 
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shifted sound wave to destructively interfere with the source wave [Lueg, 1937]. Various 
other attempts have been made using similar techniques, with those involving multiple 
cancellation sources enjoying the best results. 

The availability of faster and better processors resulted in the use of digital signal 
processing techniques to actively control sound, and for the first time in the 1970s, 
structural vibration. By the 1980s, there were several designs of active vibration control 
for systems such as the classical spring-mass damper system [Burdess, 1985] vibration 
isolation platforms [Takagami, 1985], flexible beams [Schafer, 1985], and others. The 
new techniques depended heavily on the ability to model the physical structure, sensors 
and actuators involved. This approach suffered from the inability to accurately model the 
system, or from changes in the system or environment after the model was determined. 
Many physical systems are too complex to model accurately using analytical methods, 
leading to simplifications that could severely impact the effectiveness of the control 
system. One possible answer to this problem was to use an adaptive control system, one 
that tunes itself to the system and adapts itself to changes in the system or environment 
automatically. 

Control of optical beams has only relatively recently been under consideration in 
the research community. As the amount of bandwidth necessary to transmit information 
has increased exponentially over the last decade, the use of laser beams for 
communication has become increasingly attractive. Directed energy has always been the 
plaything of science fiction, but now with the advent of new powerful lasers, the fantasy 
could become reality. However, unlike a radio wave, which spreads in a spherical 
pattern, or a missile, which can make in-flight corrections, the laser beam only goes in the 
direction you point it. Even if you can point it accurately enough, you have to deal with 
the vibrations that are inherent in mechanical systems. Jitter, as was stated earlier, 
reduces the intensity of the beam at the target, whether the beam is used as a 
communications system or a weapon. 


3 



2. Classical Control 

A straight-forward method of controlling jitter is to use the classical means of 
control such as Proportional-Integral-Derivative (PID) control or a Linear-Quadratic- 
Regulator (LQR). The PID controller works well for low frequency (less than 10 Hz) 
disturbances on the LJC, but becomes less capable at the higher frequencies. An LQR 
controller works well for controlling even at high frequencies (less than about 200 Hz on 
the LJC), but requires full state feedback (not normally available) or a means to estimate 
the states not observed [Ogata, 2002, pp. 897-910]. The LQR controller requires an 
accurate model in order to solve for the gains necessary to control the system. The ability 
to effectively control the system depends on the analytical or numerical model and the 
estimation filter. As these may be inaccurate at frequencies not allowed for, or changes 
in the environment that are not expected, the LQR may require “tuning” in order to 
maintain the desired accuracy. This tuning would not be desired in a weapons system or 
a space borne autonomous platform. 

McEver and Clark used a Linear Quadratic Gaussian (LQG) design to actively 
control jitter in an experiment designed around Fast Steering Mirrors, accelerometers and 
microphones. The disturbance was acoustically induced by loud speakers and control 
was attempted using feedforward from a microphone or the accelerometers to the FSM as 
well as feedback to the FSM using a Position Sensing Detector (PSD) at the target. 
Acoustic means of controlling the jitter were also attempted using a secondary 
loudspeaker and feedback from a microphone located near the control mirror. The FSM 
employed by McEver used a piezo-electric actuator, unlike the FJC FSMs which use 
voice coils to move the mirror. McEver’s experiment showed that using feedforward 
FSM control with the microphone resulted in a reduction from 4.4 to 2.5 rms fd rad 
broadband for the y-axis of the PSD. He did not specify any reduction in the x-axis. 
Feedback FSM control using the PSD reduced the x-axis jitter from 2.51 to 1.38 rms 
H rad and the y-axis jitter from 7.0 to 2.7 rms fu rad. The use of Acoustic means to 
control the jitter did not result in any reduction [McEver, 2001]. 
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Glaise, Anderson and Janzen conducted similar experiments as McEver, et al 
using acoustically induced jitter and LQG control. Their FSM also used piezo-electric 
actuators. Feedback control was attempted using a secondary loudspeaker as an actuator 
as well as the FSM. The feedback was provided by a PSD. The use of a secondary 
loudspeaker did not result in appreciable improvement of the jitter recorded by the PSD. 
With the FSM as the actuator, the reduction in jitter was considered good, with an almost 
complete elimination of a 30 Hz response on the PSD. The LQG controller had been 
targeted for this frequency based on the open loop response observed as a result of the 
acoustic disturbance [Glaese, 2000]. 

3. Adaptive Control 

In the early 1950s there was extensive research into methods of controlling high- 
performance aircraft using adaptive means. These aircraft operated over a wide range of 
altitude and speed. In one regime, a linear, constant-gain classical feedback system might 
work, but it would not work over the complete range of flight parameters. A more 
sophisticated controller was needed and a system of “gain scheduling” was devised that 
worked quite well. The controller thus “adapted” to its environment, allowing one 
system to be used to control a widely varying flight regime. This was a very rudimentary 
form of adaptation, but due to the limitations of computer processors at that time, was the 
best that could be done. The introduction of State-Space methods, Stability Theory and 
Stochastic Control Theory were major developments in the 1960s that aided resurgence 
into adaptive control research in the 1970s. In the late 1970s and early 1980s, proofs for 
stability of adaptive systems appeared. As processor speed and memory storage 
capability increased, the opportunity to use adaptive means improved. As a result, 
adaptive controllers started to appear commercially in the early 1980s [Astrom, 1995]. 

Adaptive control has also been used in the control of optical beam jitter. 

Skormin, Tascillo, and Busch developed a computer simulation in 1995 of an airborne to 

satellite optical link in which the use of a self-tuning regulator (STR) as well as a 

Filtered-X Least Means Squares (FXLMS) controller was used to mitigate the effects of 

jitter on the optical beam. The simulation shows that adaptive feedforward vibration 

compensation can be used to minimize induced jitter [Skormin, 1995]. In 1997, Skonnin 
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and Busch proposed the use of Model Reference Control for jitter reduction. 
Experiments were conducted using a specially designed high bandwidth FSM and a 
commercially available low bandwidth FSM. In each case, significant reduction (on the 
order of 20 dB) in acoustically induced jitter was achieved [Skormin, 1997a], Also in 
1997, Skormin, Tascillo and Busch demonstrated the use of a self-tuning regulator in an 
acoustically induced jitter rejection experiment using FSMs (voice-coil actuated), 
accelerometers and PSDs. The data showed that the STR was superior to classical 
feedback control in frequency ranges above about 300 Hz. [Skormin, 1997b]. 

a. Least-Mean-Square (LMS) Method 

Widrow and Hoff proposed the FMS algorithm in 1960 [Widrow, I960]. 
The FMS method uses a Finite Impulse Response (FIR) filter to generate a control signal 
from a reference input. The FIR gains are adjusted by the algorithm to minimize the 
Mean Square Error (MSE) of the output. The algorithm is elegant in its simplicity, 
requiring neither measurements of the pertinent correlation functions, or matrix inversion. 
The FMS algorithm was originally designed for signal processing, and is the standard 
against which other linear adaptive methods are measured [Haykin , 2002, p. 231]. When 
used as an adaptive noise or vibration canceller, the effects of the secondary path (the 
actuator) must be taken into account to improve perfonnance. This led to the Filtered-X 
version (FXFMS) of the algorithm in which the reference signal is filtered by a model of 
the secondary path. Widrow et al derived the algorithm in 1981 [Widrow, 1981] in the 
context of adaptive control and Burgess independently derived it for use in noise control 
[Burgess, 1981]. The FMS algorithm uses a reference signal for removal of the 
unwanted disturbance. This reference signal has been the subject of many papers. One 
example that is applicable to the rejection of colored noise in a laser beam is the use of 
amplitude optimization to speed up convergence for the multiple frequency case of an 
ANC system [Kuo, 1998b]. . In this paper, Kuo et al, discuss the use of phase 

compensation to reduce out of band overshoot for a single frequency ANC system. As 
we shall show, the phase, amplitude and bias of the reference signal are critical to 
convergence in a targeting scenario. 
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The FXLMS algorithm has also been used in the control of optical beams. 
As mentioned earlier, Skonnin, et al, used it in the simulation of an airborne platform to 
satellite link. Glaese, et al, used it in the experiment referenced earlier in which LQG 
control was used. In Glaese’s experiment, FXLMS was implemented with the FSM as an 
actuator. The performance varied, depending on the reference source used, from a 32 to 
47 dB reduction in jitter of an 80 Hz tone induced by the loudspeaker. The broadband 
reduction varied from 7.4 to 9.5 dB [Glaese, 2000]. Boelitz, et al conducted airborne 
experiments using the FXLMS algorithm to reject disturbances from the vibration of the 
airframe and mounted optical bench during actual flight conditions. These experiments 
resulted in varying degrees of performance, from a 25% to 99% reduction in the 
Normalized Correlated Power of the track error. In some cases the jitter actually 
increased. However, due to the complexity of the experiment, data is still being analyzed 
and the increase in jitter could be due to a loss of fine track mode during the data run 
[Boelitz, 2003] 

b. Adaptive Lattices 

In 1978, Griffiths, (1978) and in a separate paper, Makhoul, (1978) 
proposed the use of a Lattice Filter to solve a joint-process estimation problem that is 
optimal in the mean-square error sense. The Lattice Filter develops a vector of backward 
prediction errors which are orthogonal to each other, simplifying the solution to the 
problem significantly, resulting in parameters that converge faster than the normal linear 
transversal filter used in the LMS algorithm [Haykin, 2002, p. 543]. The filter developed 
by Griffiths is a member of the stochastic-gradient family of adaptive filters, of which the 
LMS filter is a part. These lattice filters are, however, computationally extensive, and in 
some cases, (specifically filters of higher orders), prohibitively so [Watkins, 2004]. Jiang 
and Gibson developed a multi-channel, unwindowed lattice filter in which the recursive 
least-squares (RLS) estimation method is used. The RLS method is more exact than the 
stochastic gradient method, but requires the inversion of the error-covariance matrix. 
Jiang and Gibson’s lattice filter uses a modified Gram-Schmidt procedure to eliminate the 
necessity of inverting any matrices [Jiang, 1995]. Boelitz, et al used Jiang and Gibson’s 
lattice filter in an airborne experiment called ABLE ACE to control structural vibrations 
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on an optical train. The results show a reduction in the Normalized Correlated Power of 
the track error from 16% to 95% [Boelitz, 2003], 

C. THESIS OVERVIEW 

Chapter II presents an outline of the basic theory used in the control algorithms 
that are implemented in the course of this research. Since the Linear Quadratic Regulator 
is widely known and used, only the basics of this algorithm are discussed. The LMS 
filter and its application to control theory will be reviewed, as well as the development of 
the Filtered-X modification. Additionally, the effects of a phase-shift of the reference 
signal are explored, as well as the effects of changing the amplitude of the reference 
signal. The Gradient Adaptive Lattice is then derived and its implementation considered. 

Chapter III will discuss the experimental setup of the Laser Jitter Control testbed. 
The sensors and actuators, as well as the computer control system and support structure 
will be explained. Chapter IV details the system identification experiments and results. 
A state-space model of the system will be derived, and the system used to calibrate the 
LJC will be detailed. 

In Chapter V, an analogous model will be derived using a simple mass-spring- 
damper system. This model will be used for investigations into the characteristics of the 
different control schemes. The analytical solution to the model will be presented and 
numerical simulations analyzed to determine the effectiveness of the different classical 
and adaptive control methods. 

Chapter VI will discuss the disturbance rejection experiments attempted on the 
LJC. The experiments are categorized into three groups; those involving an external 
disturbance to the system, those in which the platfonn containing the control system is 
disturbed, and those in which both external and platfonn disturbances occur. Chapter VII 
will provide conclusions, and recommendation for future work. 


8 



II. REVIEW OF THEORY 


A. LINEAR QUADRATIC REGULATORS 

The LQR is first developed to investigate how classical control algorithms handle 
broadband and narrowband disturbances. The system to be controlled must be modeled 
mathematically to determine the optimal gains. The LQR requires full state feedback, 
which, if not available, must be estimated. A Kalman filter or other observer may be 
used to estimate the states not measured. The linear-quadratic optimal gains (K) are 
calculated to minimize the following cost function: 


J = J (x T Qx + u T Ru ) dt 

o 


( 2 . 1 ) 


The control law is: 


u = -Kx (2.2) 

In order to obtain the optimum gains, a state-space model of the system must be 
determined. The state-space model is of the form: 


x = Ax + Bu 
y = Cx + Du 


(2.3) 


where x is the state vector and the matrices A and B describe the system, y is the output 
of the system, and the matrices C and D describe the effect of the states and the input on 
the output. In equation (2.1), the matrices R and Q are weighting matrices, to provide 
information on the relative importance of each state and the input to the control effort. 
Therefore, in order to optimally detennine the gains for a specific response, one must 
understand the relationship between the state vector and the physical quantities involved. 
Accordingly, the following model was developed for the LJC. 

A simple second order transfer function of the Newport mirror about one axis is 
defined: 
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(2.4) 


H m {s) 


(o : 


+ 2Co) n s + o/ n 


where (O n and £ are experimentally detennined. A first order system was added to 
model the optical sensor system: 


H d {s) = 


1 

Ts +1 


(2.5) 


where T is the time constant for the PSD system. The value for T was determined from 
the data given for the optical sensors. The resulting State-Space set of equations for one 
axis is given: 
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In equation (2.6), V p is the voltage at the optical sensor, Gp is the optical sensor 
gain, V/m, and D m is the distance from the control mirror to the sensor, m. 0 is the 
rotation about the control mirror’s x axis, Gm is the control mirror gain, rad/V, and V,„ is 
the voltage input to the control mirror. D/, is the distance from the disturbance source to 
the detector and 0 d is the rotation about the laser beam’s x or y axis at the disturbance 

source. 6 d is a pseudo-random variable, considered to be band-limited white noise plus 

any narrow-band disturbances. The state-space system defined above is used to solve for 
the optimal gains for the LQR, with increased weighting in the matrix Q in equation (2.1) 
on the voltage at the detector, V p . 

B. THE LEAST MEANS SQUARES (LMS) ALGORITHM 

The LMS algorithm is one of the simplest yet robust adaptive algorithms. In the 
LMS algorithm, we adjust the tap gains based on the response of the system to (1) the 
error, (2) a reference signal correlated with the disturbance, and (3) the control input. The 
algorithm uses the method of Least Squares to find the optimum values for the tap gains. 
In particular, the algorithm relies on predicting its next input, which is simply the 
disturbance in the case of laser jitter control, to optimize the tap gains. The error - the 
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difference between the predicted signal and the system output - is then used to recalculate 
the gains that minimize the error in return. For the LJC Testbed, the feed-back or target 
detector is used to provide the error signal, and the feed-forward detector provides the 
correlated disturbance input signal. This type of control algorithm not only calculates the 
necessary gains, but also identifies the system, simplifying the requirement to 
mathematically model the system. The type of adaptive control employed in the 
experimental setup utilizes predictors. One type of forward predictor is the transversal or 
ladder filter, as shown in Figure 1. One may also conduct backward prediction, i.e. 
calculate what an input was in the past, given the n current and past inputs. 



Figure 1 Transversal Filter 

The derivation that follows is modeled after Kuo (1996), modified for the optical 
jitter control system. Define an input vector at time n as: 

x(n) = [x(n) x(n- 1)--- x(n-M + \)] T (2.7) 

where M is the order of the filter, (number of weights), and the superscript T denotes 
transpose. Furthermore, the over-bar designates a vector. Define the weight vector at 
time n as: 

Mn) = K ( n ) w l (n) ■ ■ ■ w M _ x (n)f 
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( 2 . 8 ) 






The output y(n) of the transversal filter in Figure 1 is then: 


v{n) 



(n)x(n) 


(2.9) 


This output is compared with a desired signal, d(n) which results in the error equation: 


e(n) = d(n)-y(n) 

= d(n)-w (n)x(n) 


( 2 . 10 ) 


Our objective will be to minimize the Mean Square Error, £ ( n ) by adjusting the tap 
gains, w(n). We will restrict the statistics of x(n) and d(n) to be stationary and the 
vector w(n) to be a deterministic sequence. Define the Mean Square Error, (MSE) by: 


%(n) = E e (n) 


( 2 . 11 ) 


where E[ ] denotes the expected value. Substituting equation (2.10) into (2.1 l)gives: 

£(n) = E^d 2 (n)~^-2E d(n)x(n) w(n) + w ( n)E x(n)x(n) T w(n) (2.12) 


Now define the desired-to-input cross-correlation vector and the input autocorrelation 
matrix as: 


p = E d(n)x(n) 


R = E 


x 


(n)x(n) 7 


(2.13) 

(2.14) 


where the boldface letter indicates a matrix. Equation (2.12) may now be written as: 

£(n) = E d(nY -2 p w(fi)tw(n) r Rw(/i) (2.15) 


We may now determine the optimum tap gains that minimize the MSE by 
performing a vector differentiation of £(n). This results in the Wiener-Hopf equations, 
formulated for the discrete-time filter: 
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( 2 . 16 ) 


w = R 1 p 

the superscript o denoting the optimum values for the weight vector. In order to obtain 
the optimum, we must have the values for R and p , which in practice, in real time, are 
difficult if not impossible to evaluate accurately. We must have a means to estimate R 
and p. 

It can be shown that the MSE is a quadratic function in the weight vector, and that 
— 0 

any departure from w will result in an increase in the MSE. The MSE surface has only 
one global minimum point, and at that point, the tangents to the surface must be zero. 
We may therefore use gradient-based algorithms to search out the optimum weight vector 
[Kuo, 1996, pp. 22-23]. The LMS algorithm makes use of the steepest descent method to 
find the optimum weight vector. The steepest descent method takes “steps” in the 
direction of the negative of the gradient of the surface, here the MSE surface. We update 
the weight vector as follows: 

w(n + \) = w{n)-^W£{n) (2.17) 

where ju is the step size that controls stability of the filter as well as the rate of descent 
along V^(«). The factor 2 is used to simplify the final result. From equation(2.15) we 
may obtain V £ ( n ): 

V^(n) = -2p + 2Rw(n) (2.18) 

However, as stated earlier, we need to find a means to estimate p and R. Widrow 
(1970) used the instantaneous squared error as the estimate for the MSE: 

I {n) = e 2 (n) (2.19) 

The estimate of the gradient of the MSE is thus: 

VJ(/i) = 2[Ve(/i)]e(/i) (2.20) 
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and, noting that Ve(n) = -x(n) (see equation (2.10)), we find that the update equation 
for the weight vector to be: 

w (n + i) = w(n) + jUx(n)e(n) (2.21) 

which is kn own as the LMS algorithm. The block diagram of the algorithm in a form for 
implementation on the LJC is shown in Figure 2: 



P denotes the primary plant transfer function which the disturbance must pass 
through before the output, and S denotes the secondary plant or actuator transfer function 
the control signal must pass through before the summing junction. P includes the transfer 
function of the structure of the bench and mirrors for the case of vibration of the bench, 
the transfer function of the sensors, as well as the gain effect caused by the distance the 
light beam travels through the system to the target. S must include the effect of vibration 
to the control mirror, the transfer function of the mirror and sensors, and the delays 
inherent in the digital signal processing and the computation of the control signal. In 
Figure 2 it may be seen that the output of the FIR filter does not go directly to the 
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summing junction for the error. Therefore the conventional LMS algorithm must be 
modified to ensure convergence. Taking the Z-transform of the signals in Figure 2, we 
find that the error signal is: 

E(z) = [P(z)-S(z)W(z)]x(z) (2.22) 

Thus W(z) must model the plant as: 

, \ P(z) 

=-H (2.23) 

S(z) 

in order to drive E(z) to zero. We note here that it is impossible to compensate for the 
delay caused by S if the primary plant P does not contain a delay of at least equal length. 
For the control of optical jitter, the delay in S can result in severe difficulties, due to the 
velocity of light, which is nearly instantaneous when compared with the speed of the Fast 
Steering Mirror (FSM) and control circuitry. And in general, if S is zero for some 
frequency, the controller will be unstable. Therefore, the effects of the secondary plant 
must be taken into account. 

A simple solution to the problem is to provide an inverse filter (1 /S(z)) in series 
with S(z). However, the inverse does not always exist. A second solution would be to 
place a copy of the secondary plant transfer function in the path to the updating algorithm 
for the weight vector. This is known as the Filtered-X LMS (FXLMS) algorithm and was 
derived by Widrow (1981) and Burgess (1981). The block diagram is provided below: 
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Figure 3 FXLMS Block Diagram 
The update equation for the weight vector (equation(2.21)) now becomes: 


w (n + 1) = w(n) + jux'(n)e(n) 


(2.24) 


with x'(n) defined as: 


x'(n) = s(n)* x(n) 


(2.25) 


and s being the estimate of the impulse response function of the secondary plant, the * 
denoting convolution. 

The value of the step size // has been the subject of no little study, fu effects 
both the rate of convergence as well as the stability of the algorithm. The correct 
selection of a step size is usually done with knowledge of the maximum value for 
stability, and then by trial and error. There have been a few studies done on the optimum 
value of [A for different cases, (see for example Haykin (2002), pp. 324-325) but these 
are in general hard to implement in real time. For the FXLMS algorithm, Elliott and 
Nelson (1993) have shown that the maximum step size is approximately: 


/Anax 


(«) (M + A) 


(2.26) 
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where A is the delay in the secondary path, and M is the filter length. Another means for 
determining the step size is to use the normalized LMS algorithm. In this scheme, the 
step size is varied according to the input power of the reference (or filtered reference) 
signal, /u becomes a function of the iteration, n, and may be calculated as [Kuo, 1996, 
pp. 33-35]: 

//(») = - r ( 2 . 21 ) 

x (n)x(n) 

where the parameter a is bounded by 0 < a < 2 . The normalized FXLMS (nFXLMS) 
algorithm is used in most of this research. 

Another factor in the selection of the step size is the Misadjustment, M,. The 

misadjustment is the ratio of the excess MSE to the minimum MSE. The excess MSE is 

—° 

caused by the algorithms inability to converge to the optimum weight vector, w . Since 
we are using the steepest descent method, there will be a certain amount of gradient noise 
due to the fact that the gradient of the MSE must be estimated. This noise results a MSE 
greater than the MSE that results from the Wiener-Hopf solution. M, can be shown to be 

directly proportional to /u as: 

u M 

M l P 28 ) 

^ k=l 


where A k are the eigenvalues of the autocorrelation matrix R and M is the filter length. 
As // is increased, the misadjustment increases and the excess MSE becomes greater. A 
second factor in the selection of // is the time constant, or rate of convergence. The 
average time constant (the time it takes for the MSE to decrease by a factor of 1/e from 
its starting value) is calculated as: 


1 

H 


(2.29) 
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where A av is the average of the eigenvalues of R [Haykin, 2002, pp. 271-272], From 
equations (2.28) and (2.29) we see that as n is increased, the excess MSE increases, but 
the average time constant decreases. There is a tradeoff in the selection of the step size. 

1. Effect of Phase Shifting the Reference Signal 

Let us define an input vector as: 

v(«) = [v(«) v(n- 1) ••• v(«-M + l)] (2.30) 

This vector is the actual disturbance before the measurement that results in x{n). We 

restrict v(n) at this time to be the vector of discrete samples that result from a periodic, 
zero-mean disturbance. We now change our FXLMS block diagram to include an 
Internally Generated Reference Signal (IGRS) that generates x{n). In other words, 

sensors detect v(n) and generate x(n ). Now, as discussed earlier, we note that any 

error in the optical bench’s position or the laser beams position is instantaneously 
(relative to the FSM’s velocity) measured at the optical error sensor. This effect results 
in a violation of causality. Flowever, as long as the beams position changes slowly 
enough, one can correct the error before it becomes too large. The block diagram of the 
system looks like this: 



Figure 4 


FXLMS with IGRS 
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Let v( n ) be A sin {'lnfnT : ). Then, d(n) will be PA sin [2nfnT s +(l> ;t ) where T s is 
the sample time, / is the frequency in cycles/sec and (f) p is the phase delay through the 
primary plant. Now, let the generated reference signal, x{n) be: 

x(n) = Bsin(2KfnT s +(j> G ) (2.31) 

where (p G is to be determined. The output of the FIR filter generated by the LMS 
algorithm is: 

u(n ) = IVB sin (2k fnT s +(/> G + (j) w ) (2.32) 

where (j) w is the phase delay through the FIR filter. Let the output of the secondary plant 
(the FSM and DSP circuitry) be s(n). With input u(n), s(n) will be: 

s(n) = SIVB sin (2k fnT s + <J) G + (j) w + <p s ) (2.33) 


where (j) s is the phase delay through the secondary plant. In order for e(n) to be zero, 
5 ( n ) = d (n) which results in: 


and: 


SWB = PA 


(2.34) 




(2.35) 


The function of the LMS algorithm it to adapt W and (f) w in equations (2.34) and (2.35) to 
accomplish compete rejection of the disturbance signal, i.e. 

PA 

W = — (2.36) 

SB 


and: 


0S 


(2.37) 
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Therefore, if one could determine appropriate values for W and (f) w such that the 
algorithm started the recursion with tap gains near the optimum values, the time it would 
take for the algorithm to reach steady state would be reduced. 

Now, returning to the weight vector update equation (equation (2.24)) we see that 
the rate of update of the vectors is dependent on three factors for the nFXLMS 
algorithm: 1) the value of /J, 2) the input power of the reference signal, and 3) the 

correlation between x(n) and e(n). We note that x'(n)e(n) in equation (2.24) is a 
rough estimate of the desired MSE gradient, however: 

E(vf) = E(e(n)x'(n)) = Rw(n)-p = V£ 2 (2.38) 


and it is therefore an unbiased estimate [Shan, 1988]. If we maintain // and the 
reference signal power constant, then the adaptation rate is strictly dependent on the 
correlation between x'(n) and e(n). When the correlation is the greatest, the adaptation 

rate of the weight vector is greatest, and when x'(n) is uncorrelated with e(n), the 

adaptation of the weight vector stops. Since the change in the weight vector is always in 
the direction of the estimate of the negative gradient of the MSE surface, we conclude 
that by adjusting (f) G to better correlate x '(/?) with e(n) the rate of convergence to the 
optimum weight vector may be increased. This effect decreases with increasing //, but 
the use of a large fl increases the misadjustment. We can therefore decrease the time 
constant of the process without increasing the misadjustment. 
a. The Adaptive Delay Filter 

The question now becomes what to set for a value of (/) G 1 To start with, 
let us evaluate the phase of the error: 


e(n) = d(n)-w (n)x'(n) 

= PAsin(27rfnT s +(f> p )-SiVBsin (ijlfnT + (f) G + (f) w + </> s ) 


(2.39) 
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Through the use of trigonometric identities, one may determine the phase shift of the 
error (</> e ) to be: 

. ... PAsm(</>)-SWBsm{(/> G +</> lv +(p s ) 

tan(<z>J =- 7 —r- (2.40) 

PA cos y(j ) p ) - SWB cos (<j) G + (f) w + (/) s ) 

The evaluation of (j) w is difficult, since it is being adapted by the LMS algorithm, but at 
the outset of the adaptation process, it is easy to see that (j) w =(j) G . If we set the initial 
values of w(n) such that they are all equal (nominally zero), the w{n) will just be an 
attenuated x{n) after the first adaptation step, and therefore, (j) w = (j ) G . The correct (f) G 
will be that which equals (j ) e , resulting in the optimum correlation between e(n) and 
x(n ). We may thus plot the value of tan (<j) G ) verses the tan( <j ) e ) as a function of (j) G to 
find the value that maximizes the first adaptation step, assuming we know (j) p and (j) s . 
However, at the next adaptation step, <j) w begins to change, and the correct value of (j) G is 
no longer known. We also may not know the precise value of (j) p and </> s , or their values 
may change with time due to changes in input frequency, etc. (j) G should be adaptable in 

order to take advantage of this effect. Kuo and Chung [Kuo, 1998b] developed an 
Adaptive Delay Filter (ADF) that used one variable delay element and one variable gain 
element to estimate the time delay and gain difference between the primary and 
secondary paths. However, their filter used the time delay element and gain element to 
replace the FIR filter W(z). We make use of this adaptive delay element in a different 
scheme to vary the generated reference signal phase as can be seen in the block diagram 
below: 
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Figure 5 nFXLMS with ADF 


In the weight vector update equation for the ADF, a symmetric difference 
is used to approximate the gradient to reduce the influence of bias from gradient 
estimation [Etter, 1981], In the ADF algorithm, there is only one output, the delay that 
corresponds to <f> G , A G . The update equation for A G is: 

A G [n + 1) = A G («) + /i A e(n)[x(n- A G (n)- A s -\)-x[n — A G (n)- A s + l)J (2.41) 

where jU A is the adaptation rate for the ADF , e ( n ) is the error, and A v is the delay 
through the FSM and control circuitry. ju A is bounded by [Etter, 1981]: 

0 <jU A <— (2.42) 

A lOE 

where P x is the power in the input signal. The ADF algorithm finds its main use when 
the adaptation rate for the LMS algorithm is required to be small in order to minimize the 
misadjustment. In the LJC experiment, due to the causality issue mentioned above, the 
experimental values of (j) G results in a Lead angle for frequencies above about 10 Hz. In 
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practice, A G is started out with a delay of about one period. In this manner, the ADF 
algorithm can then adapt to a “lead” angle, thus correcting the problem of (j) p being less 
than </> s . 

2. Effect of Bias on the LMS Filter 

In most laser targeting schemes, some sort of compensator is used to correct the 
bias error at the target. A second controller is then used to remove the “noise” in the 
beam. It has been noted by Widrow that an adaptable bias weight may be used to 
counteract the effect of “plant drift” in an adaptive inverse modeling situation [Widrow, 
1985, pp. 283-284]. This single adjustable weight adapts to remove any bias in the 
output of the plant. It was noted during testing on the LJC, that by adding bias to the 
reference signal, the LMS fdter adapted more quickly to the disturbance, and that the 
signal was zeroed on the target, without the use of any compensator. This effect may be 
analyzed as follows. Referring to Figure 4, let v{n) = C + A sin [lKfnT s ) where C is a 
bias error. Then the desired signal (the signal that must be cancelled) is: 

d (n) = PC + PA sin (2 nfnT s + </> p ) (2.43) 

where P is the gain through the optical bench and mirrors and (j) p is the delay, as before. 
Now, considering only the bias effect, let the generated reference signal be: 

x(n) = C'+ B sm(2nfnT s ) (2.44) 

where C’ is an adjustable bias. The output of the FSM will be: 
e(n ) = d (n)-s(n) 

r , \ 1 r (2.45) 

= PC + PA sin (In fn T s + </> p )\- [SWC ’+ SWB sin (2 xfnT s + (j) w + </> s )] 

Rearranging to place the DC components together we find: 

e{n) = [PC-SWC'] + \PAsm(27rfnT s +^ p )-SWBs,m{27rfnT s +(j) w +(j) s ) \ (2.46) 

It is seen here, that the proper selection of C’ will result in the cancellation of the DC 

component of the error. It is also noted, that if C’ is zero, the LMS algorithm will be 
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unable to completely cancel the resulting error. As with the ADF, C’ should be 
adaptable, since C may change and since W is adapting during the process to some quasi¬ 
steady state value. 

a. The Adaptive Bias Filter 

A similar process to the ADF is applied in the development of the 
Adaptive Bias Filter, or ABF. Using a one or two weight nLMS filter, the bias in the 
reference signal is adjusted to remove the DC component of the error signal. An estimate 
of C’ is used as the reference signal to the ABF. The error signal in this case is the mean 
error, which stops the adaptation once the signal is centered on the target. 



Figure 6 nFXLMS controller with ADF and ABF 

C. THE GRADIENT ADAPTIVE LATTICE 

The Gradient Adaptive Lattice (GAL) algorithm uses both forward and backward 
prediction to develop the estimate of the desired cancellation signal, which again drives 
the control mirror to cancel the effect of disturbance input. The construction of the GAL 
filter is different than the linear transversal filter in that it uses reflection coefficients vice 
tap gains due to the nature of the lattice-like structure as shown in the figure below. 
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The forward prediction error for the m th stage is given the symbol f m ( n). The 
backward prediction error, b m is given likewise. K is the reflection coefficient for the m th 
stage and h m provides the gains to construct the desired response, y m . The h m are 
calculated using the method of least squares, with the b m as the regressor and the h m as 
the parameters to be solved. The lattice filter results in a highly efficient structure in that 
it is modular. Each stage only requires the output of the previous stage, and an increase 
in the order of the filter only requires adding additional stages. Jiang, et al [Jiang, 1995] 
and Chen [Chen, 2001] discuss active noise cancellation techniques using algorithms that 
employ the lattice filter. The same concept was used in this study, but the algorithm from 
Haykin [Haykin, 2002, pp. 536-543] was used due to its simpler implementation. 
Although modular in its implementation, the GAL is computationally expensive, as the 
calculation of the reflection coefficients is very complex. 

In brief, the lattice filter reflection coefficients are detennined by minimizing the 
cost function: 


J = -E 
2 


fm 00 ‘ + K 00 


(2.47) 
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Using the input-output relationships for the f n and b m one can determine the 
optimum K by differentiating the cost function with respect to K . This leads to the Burg 
estimate [Burg, 1968] for the values of the optimum reflection coefficient: 

22x, 

= -^7^---- (2 48) 

z i.Ci(*)r+h-,(*-i)r 

k =1 V 

where the * denotes conjugation. The determination of the h m for the desired response 
estimator uses the normalized LMS algorithm, with the input being the backward 
prediction error, b m , as the reference signal. It can be shown that the b m are all orthogonal 
to one another [Haykin, 2002, pp 172-173]. The advantage to the GAL is that since the b m 
are all orthogonal, the input correlation matrix ( R ) is a diagonal matrix, making the 

calculation of the h m relatively simple, since by the Weiner-Hopf equations, h = R p . 
We can thus directly calculate the optimum tap weights for the desired response 
estimator. This allows for a rapid convergence of the algorithm to the optimum output. 

Another feature of the GAL is its use of a time varying step size in the 
determination of the reflection coefficients. The step size is inversely proportional to the 
energy in the prediction error. A small value for the prediction error indicates that the 
predictor is providing an accurate model. This normally works well in that when the 
prediction errors are small, the algorithm uses a large step size and the filter can respond 
rapidly to changes [Haykin, 2002, p. 542]. However, if the input has a random 
component, as may be expected with a satellite laser relay system, the prediction errors 
become large and the algorithm responds much more slowly, which may be a detriment 
to the GAL in its use with the LJC in controlling “colored noise.” 


26 



III. EXPERIMENTAL SETUP 


A. LASER JITTER CONTROL TESTBED 



The Laser Jitter Control (LJC) Testbed is located in the Spacecraft Research and 
Design Center - Optical Relay Mirror Lab at the Naval Postgraduate School in Monterey 
CA. The components are mounted on a Newport Optical Bench which can be floated in 
order to isolate the components from external vibrations. The idea is to simulate a 
satellite or vehicle “bus” which houses an optical relay system. The laser beam 
originates from a source and passes through a Disturbance injection Fast Steering Mirror 
(DFSM). The DFSM corrupts the beam using random or periodic disturbances 
simulating the disturbances that might occur as a result of the transmitting station or the 
tip and tilt the beam may suffer as it passes through the atmosphere. A vibration isolation 
platform is used to mount the relay system, and to isolate the relay system from the 
optical bench. The relay system platform may be vibrationally disturbed by a 5 lb inertial 
actuator, simulating on-board running equipment such as control moment gyros, reaction 
wheels, cryo-coolers, etc. The incoming beam is split and reflected off the platfonn to a 
Position Sensing Detector (PSD) in order to obtain a reference signal that indicates the 
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on-board and injected disturbances. It is recognized that this is not a normal means to 
measure the disturbances as one would not be able to mount a PSD off the ‘bus’. 
However, this reference PSD may be seen as simulating an on-board Inertial 
Measurement Unit (IMU), which is normally available in satellites with an optical 
payload. Using an accurate and sensitive IMU and a PSD, one could detennine a similar 
reference signal. This setup allows a comparable measurement without the added cost of 
an IMU. A control FSM, designated the Receive FSM or RFSM, is used to correct the 
disturbed beam. The corrected beam is then reflected off the platfonn by a second 
folding mirror to the target PSD. 



Figure 9 Laser Jitter Control Testbed 


B. FAST STEERING MIRRORS 

The Fast Steering Mirrors are the heart of the LJC. They are used to rapidly and 
accurately direct the laser beam through the system. The FSMs in the LJC use voice 
coils to position the mirrors in response to commanded inputs. The LJC uses two 
different FSMs, one by the Newport Corporation, and one by Baker Adaptive Optics. 
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1. Newport Fast Steering Mirror, Model FSM200 



Figure 10 Newport Fast Steering Mirror 


The Newport FSM is used as the control mirror (RFSM) for all experiments 
conducted during this research. The mirror comes with its own controller, the FSM- 
CD100, capable of both open loop and closed loop operation. The controller also 
provides an output of the mirror’s angular position about each of the axis. In these 
experiments, the controller is configured in the open loop mode, with control inputs 
provided from the computer control system. The FSM200 is equipped with a two inch 
diameter mirror. The mirror substrate is mounted to its support by means of a flexure 
suspension, which does not require bearing surfaces, eliminating stiction and wear. The 
suspension for the FSM200 has a very small spring constant to minimize the current (and 
the subsequent heating) necessary to move the mirror. As a result, the mirror “hangs” 
when in the powered-off condition, and must be powered on before aligning the laser 
beam. It also must be powered on before shaking the platfonn to prevent damage to the 
suspension system. The actuator voice coils are mounted behind the mirror, and operate 
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in a push-pull manner, rotating the mirror about the axis that bisects each pair. The axes 
of the FSM200 are aligned with the horizontal and vertical directions as referenced to the 
mirror housing. The FSM has a control bandwidth of about 800-1000 Hz [Newport, 
2002]. The manufacturer’s specifications and operating limits are provided in Appendix 
A. 

2. Baker Fast Steering Mirror, Model “Light Force One” 



Figure 11 Baker Fast Steering Mirror 


The second FSM used in the LJC is from Baker Adaptive Optics. The Light 
Force One is a one inch diameter mirror used for disturbance injection (DFSM) in these 
experiments. The Baker mirror uses a stiffer spring to mount the mirror substrate, which 
allows the mirror to remain in position when powered off. The Baker mirror comes with 
a small driver for positioning the mirror; however there is no closed loop option and 
mirror position is not available. The control bandwidth for the Baker mirror is about 3 
kHz. Specifications are provided in Appendix A. 
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c. 


POSITION SENSING DETECTORS 



Figure 12 Position Sensing Module 


The laser beam optical position sensors, known as Position Sensing Modules or 
PSMs, were purchased from OnTrak Photonics Inc. and have a detection area of 10mm x 
10mm. Each duolateral, dual axis PSM (model PSM2-10) requires an amplifier, the OT- 
301 to output the x and y position (in volts) of the centroid of the laser beam spot. The 
combination of amplifier and detector is called a Position Sensing Detector (PSD). 
Various beam splitters and folding mirrors are used to direct the beam as desired. The 
detectors have a minimum sensitivity of 50-100 /d W, which drives the selection of the 
beam splitters and the determination of laser power. Their frequency response for the 
gain settings normally used is 15 kHz. The optimum beam diameter for accurate 
detection on the PSM is between 1 and 3 mm, and the maximum allowed intensity should 
be less than 300 W/cm 2 . The output range of the PSM2-10/OT-301 detector is ±10 V and 
the OT-301 amplifier has a noise level of 1 mV. The minimum resolution of the PSM2- 
10 with the OT 301 amplifier is therefore 0.5 fu m [Hunnicutt, 2003]. The specifications 
and drawings are provided in Appendix A. 
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D. NEWPORT VIBRATION ISOLATION PLATFORM 



Figure 13 Newport Vibration Isolation Platform 


The RFSM, beam splitters and folding mirrors are mounted on a Bench Top 
Pneumatic Newport Vibration Isolation Platform. This platform allows the control 
system actuators and optical path to be vibrated independent of the Optical Bench. The 
breadboard, which is self-leveling, rests on four pneumatic isolators. 
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The response of the platform to vibrations from the workbench surface is shown 


below: 
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As noted above, the Isolation Platfonn was designed to isolate the breadboard 
from the vibrations of the workbench surface. In our experiment, we wished to vibrate 
the breadboard, and isolate the workbench surface from the vibrations. This concept was 
discussed with the Newport technicians, and they agreed that the system would work in 
reverse, as long as the vibrations were of such small amplitude as to not cause the self¬ 
leveling system to attempt to correct. The vibrations we intended to induce were on the 
order of 10 to 500 microns, which was well outside the performance specifications of the 
leveling system, and thus it would not correct the tilt of the platfonn induced by the 
mounted inertial actuator. 

E. CSA ENGINEERING INERTIAL ACTUATOR 

In order to vibrate the platform at desired frequencies, an inertial actuator is 
mounted on the vibration isolation platfonn. This actuator is a CSA model SA-5, capable 
of delivering a rated force of 5 lbf, in a bandwidth of 20 to 1000 Hz. The SA-5 has a 
resonance at about 40 Hz. The “shaker” may be mounted in any orientation. For our 
experiments, it was mounted next to the last folding mirror before the beam exits the 
platform, in the vertical direction, to provide maximum x and y axis motion. A separate 
bracket was designed to allow mounting in the horizontal direction as well. The 
manufacturer’s specifications and drawings are provided in Appendix A. 
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F. KISTLER 3-AXIS ACCELEROMETER 

A Kistler model 8690C10 accelerometer is employed to measure the disturbances 
generated by the shaker on the platform. The accelerometer is driven by a Kistler 
Piezotron Coupler, model 5124A. 

G. COMPUTER CONTROL SYSTEM AND SOFTWARE 

The computer control system is based on MATLAB, version 6.1 release 13 with 
SIMULINK, and the xPC Targetbox, all from the Mathworks. The main computer for 
control implementation and experiment supervision is a 2.4 GHz Dell Dimension 8250 
with 1 Gbyte of RAM. The xPC Targetbox has the ability to accept 16 differential inputs 
and provide 4 analog outputs. The Targetbox is a Pentium III class computer running at 
700 MHz. The Disturbance Computer, currently controlling the shaker, uses an AMD 
Athlon processor running at 1.4 GHz, with 256 Mbyte of RAM. dSPACE ver 3.3 with 
ControlDesk ver 2.1.1 is used to interface with the shaker. Future plans envision the use 
of two control mirrors. This will require the disturbance computer to control the DFSM 
as well as the shaker, in order to free up analog outputs from the xPC Targetbox. 
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Figure 16 Signal Flow Diagram 


The Control PC is used for several functions. Coding for control algorithms is 
performed in the MATLAB/SIMULINK environment. Some specialized functions, such 
as the lattice algorithm, require SIMULINK “S-functions”. These are coded in C and 
compiled using the installed C++ compiler. xPC Target software compiles, generates 
code and downloads to the Targetbox the SIMULINK model for use in the control 
scheme. The Targetbox then runs the code automatically. The control PC provides 
supervisory functions during code execution, and allows user interface to the running 
code through the SIMULINK model. The disturbance computer uses the ControlDesk 
software to input the desired frequency, multiple frequencies, or random signals to the 
shaker through a Kepco Power supply. The output of the target PSD, OT2, is also routed 
to the disturbance PC through the dSPACE interface to provide a real-time plot of the 

beam position. The versions of software used for the LJC are listed in Appendix B. 
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IV. SYSTEM IDENTIFICATION EXPERIMENTS AND 

CALIBRATION 


A. SAMPLE RATE 

The frequencies of interest for the control of optical beam jitter vary from band- 
limited white noise at up to 1 kHz, to multiple narrow band frequencies in the range of 10 
to 200 Hz. Due to the control bandwidth of the Newport mirror (800 Hz) and some trial 
runs, we determined that the mirror could handle up to one-fourth of its control 
bandwidth and maintain a reasonable control effort (i.e. not exceed the mirror’s current 
limitations). We therefore chose to limit the random input to 200 Hz band-limited white 
noise. Since we were attempting to control a “noisy” signal, we decided to use the fastest 
sampling rate possible without overloading the CPU. Once the control algorithms and 
SIMULINK models were designed, we detennined a sample rate of 2 kHz was possible. 
It is unlikely that any noise or frequency components he above the Nyquist frequency of 
1 kHz for our system. Additionally, any filtering of the input signal will necessarily incur 
a phase delay, which in the case of random noise at 200 Hz could adversely affect the 
ability to mitigate the jitter in the optical beam. For these reasons, an anti-aliasing filter 
was not used. 

B. STEP, IMPULSE AND FREQUENCY RESPONSE OF THE FSM 

OPTICAL SENSOR SYSTEM 

As a first attempt to model the system, Bode diagrams for the X and Y axis of the 
Newport FSM and Baker FSM were determined. The Baker plots used the voltage-in as 
the input, and the output was taken from a PSD located near the mirror, as the Baker 
mirror has no position sensor. The shape of the curve, especially the response at 
frequencies above 1 kHz, was discussed with the manufacturer, who reported that the 
response was typical of the mirror as designed. The speed of response of the Baker 
mirror is such that a very short (less than 4 microseconds) “spike” in the output of the 
D/A converter of the xPC Targetbox is picked up by the Baker mirror. The spike occurs 
at the “zero-crossing” of the signal, and is believed to be a problem with the circuit board 
in the Targetbox. 
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Baker X axis Bode plot using 400 avgs. Sampling interval 0.0002 secs. Freq res 0.05 Hz. NyFreq 2500 Hz 
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Figure 17 X axis Bode Plot of the Baker Mirror 

Baker Y axis Bode plot using 100 avgs. Sampling interval 0.0002 secs. Freq res 0.2 Hz. NyFreq 2500 Hz 
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Figure 18 Y axis Bode Plot of the Baker Mirror 
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The spike currently prevents effective use of the Baker Mirror as a control mirror, 
and forces its use as the DFSM. The Newport mirror does not pick up this spike, and is 
not affected by it. This is probably due to the lower bandwidth of the Newport mirror. 


Newport X axis Frequency Response Function using 400 avgs. Sample Rate 5 kHz. Freq res 0.033333 Hz. NyFreq 2.5 kHz 
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Figure 19 X axis Bode Plot of the Newport Mirror 


39 



















































Newport Y axis Frequency Response Function using 400 avgs. Sample Rate 5 kHz. Freq res 0.033333 Hz. NyFreq 2.5 kHz 
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Figure 20 Y axis Bode Plot of the Newport Mirror 


The Newport Bode plots have the 3 dB down points marked to show the 
bandwidth of the mirror. The detector output is plotted also, indicating the response of 
the detector. The detector curve shows the effect of the additional pole by comparing the 
slope of the magnitude line for the mirror with that of the detector at frequencies above 1 
kHz. The responses for the Newport mirror are typical of that for a highly damped, 2 nd 
order system. 

Using the frequency, step, impulse and sinusoidal responses of the Newport 
mirror (RFSM) and target detector, a State-Space model was determined using classical 
analysis, similar to that given in Ogata [Ogata, 2002, pp. 224-233]. The parameters were 
adjusted to give the best fit for the sinusoidal response in the frequency range of 50 to 
100 Hz, since we are primarily concerned with the ability of the adaptive systems to 
manage the vibrations of the platfonn. The State-Space model was constructed to allow 
for cross-coupling between axes. A comparison of the model to the actual system is 
given below for the X axis of the RFSM. 
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Newport X axis Response 
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Figure 21 


Impulse Response, Newport X axis 
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Figure 22 Step Response, Newport X axis 
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Newport X axis Response 



Figure 23 50 Hz Sinusoidal Response, Newport X axis 

Newport X axis Response 



Figure 24 100 Hz Sinusoidal Response, Newport X axis 
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Newport X axis Frequency Response Function using 400 avgs. Sample Rate 2 kHz. Freq res 0.033333 Hz. NyFreq 1 kHz 




Figure 25 Frequency Response, Newport X axis 


Although the impulse and step response of the model do not compare as favorably 
with the output from the mirror as may be desired, the parameters are determined with 
emphasis on the frequency response characteristics. 

C. STATE-SPACE MODEL OF THE SYSTEM 

Using equation (2.6), and expanding it to cover both axes, the State-Space model 
for the LJC Testbed is given below. The parameters used in the model response given in 
the figures above are provided in Table 1 
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(4.2) 


Name 

Variable 

Value 

PSD Voltage, Y axis 

y w 

-10 to +10 volt 

PSD Voltage, X axis 

v * 

-10 to +10 volt 

Rotation about X axis of RFSM 

o x 

-26.2x10' 3 to +26.2x10' 3 rad 

Rotation about Y axis of RFSM 

6y 

-26.2x10' 3 to +26.2x10' 3 rad 

PSD response Time 

T 

67 xlO' 6 sec 

Geometry factor, Y axis 

gte y 

0.405 (typical) 
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Name 

Variable 

Value 

Geometry factor, X axis 

&0X 

0.535 (typical) 

PSD Calibration, Y axis 

GPr 

2xl0 3 V/m 

PSD Calibration, X axis 

G Px 

2xl0 3 V/m 

Distance from RFSM to Target 

Dm 

1.245 m 

RFSM Damping Ratio, X axis 

c 

0.90 

RFSM Damping Ratio, Y axis 

c 

0.90 

RFSM Bandwidth, X axis 


5655 rad/sec (900 Hz) 

RFSM Bandwidth, Y axis 

CO y 

5184 rad/sec (825 Hz) 

Mirror Calibration Factor 

Gm 

52.4x10' 3 rad/volt 

Voltage input to RFSM, X axis 

V 

mx 

-10 to +10 volt 

Voltage input to RFSM, Y axis 

v my 

-10 to +10 volt 

Cross-coupling Factor, X axis 

A 

-2 xlO’ 2 

Cross-coupling Factor, Y axis 

A 

-2 xlO' 2 

Table 1 List of Variables for t 

he State Space Model 


In equation (4.1), the geometry factors, gt 0 and gt 0x are used to correct for the 

45 degree angle the RFSM is mounted at with respect to the incoming beam, as well as 
the alignment of the folding mirrors and beam splitters. The geometry factors vary 
depending on the height of the vibration isolation platfonn, which varies slightly with 
pressure, but is measurable on the PSDs. The cross-coupling factor was chosen by 
comparison of the output of the opposite axis with input on one axis. The RFSM 
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bandwidth and critical damping ratio were chosen by comparison of the model output 
with the actual frequency response (Bode diagram) and the sinusoidal response in the 
frequency range of interest (50 to 100 Hz). All other numbers are from manufacturers 
data. The factor “-1” in equation (4.2) takes into account the wiring difference between 
the feedback detector (OT3) and the target (OT2) and reference signal (OT1) detectors. 
OT3 faces OT2 and OT1 and therefore its X axis is reversed in direction from OT2 and 
OT1. 

D. FREQUENCY RESPONSE OF THE PLATFORM 

As part of the investigation into controlling the vibrating platform, we need the 
frequency response of the RFSM/platfonn system. A linear chirp signal was used to 
vibrate the inertial actuator, from 10 to 500 Hz, with a period of 1 sec and amplitude of 1 
volt peak to peak. The Power Spectral Densities calculated in this research are all 
determined by the pwelch function in MATLAB which uses a Hamming Window to 
shape the data. The Y axis response on the PSD (rotation of the X axis of the mirror) is 
shown below: 


PSD Plot of Platform; Shaker at (28,17,Y), accel at (31,17). 



Figure 26 Power Spectral Density of the Platform 
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Figure 26 shows the response of an accelerometer placed next to the inertial 
actuator and the sensor response at both the target (OT2Y) and reference (OT1Y) PSDs. 
Power Spectral Density Plots were obtained using accelerometers at the top and bottom 
of both the Newport mirror and the last folding mirror on the platform, labeled FM2 and 
compared to the accelerometer located next to the inertial actuator. 


PSD Plot of Newport; Shaker at (28,17,Y), accel at (27,6) 



Figure 27 Power Spectral Density of the Newport Mirror 
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Figure 28 Power Spectral Density of Folding Mirror 2 


For the frequencies of concern, approximately 50 to 100 Hz, we note that there 
are three major areas of interest. The peak at about 40 Hz in all three figures above is 
related to the shaker resonance. The peaks at 56 and 67 Hz (Figure 26) appear to come 
from Folding Mirror 1 (FM1) and only show up on the OT1 (reference signal) response. 
The peak near 100 Hz comes from FM2 (Figure 28). A look at the Frequency Response 
Function (FRF) using the control voltage to the shaker power supply as the input, and the 
sensor data (Y axis) as the output will serve to illustrate these points: 
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Y axis Frequency Response Function using 400 avgs. Sample Rate 2 kHz. Freq res 0.033333 Hz. NyFreq 1 kHz 




Figure 29 Frequency Response Function, shaker input to OT1Y, OT2Y 

The FRF indicates that both target and reference PSDs develop a phase shift 
starting at 30 Hz. OT1Y, (reference signal Y axis) however, is about 90 degrees out of 
phase with OT2Y (target signal Y axis) between 40 and 50 Hz. OT1Y is 180 degrees out 
of phase with OT2Y between 70 and 80 Hz. As adaptive filters require a reference 
signal, these phase shifts will become important as we attempt to control the vibrations. 
It is noted, that even if one were to use the accelerometer located at the shaker as a 
reference signal, there is still a phase shift between the accelerometer and the target (30 to 
50 Hz). The table below lists the significant frequencies excited between 0 and 200 Hz, 
and the probable source for each. 
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Frequency 

Source 

38 Hz 

Inertial Actuator resonance, Folding Mirror 1 and 2, 

40 Hz 

Platform resonance 

43 Hz 

Folding Mirror 1 and 2 

55 Hz 

Folding Mirror 1 

67.5 Hz 

Folding Mirror 1 

100 Hz 

Folding Mirror 2 

196 Hz 

Platform 


Table 2 Frequency Sources 


1. Damping of the Inertial Actuator/Platform System 

The damping of the Inertial Actuator/Platform System was determined using the 
half-power points of the main resonance at 40 Hz. Figure 30 below is the Frequency 
Response of the Platform to a swept sinusoid forcing function, 35 to 47.5 Hz, to the 
inertial actuator. An HP 33120 Wavetek Signal Generator was as the source for the 
inertial actuator. An HP 35665A Digital Signal Analyzer was used to analyze the signal. 
The output of the system is the laser position as recorded on the target PSD, OT2. The 
damping calculation is shown in Figure 31 and is approximately 0.20% at 40 Hz, which 
matches the Digital Signal Analyzer determination of the damping at the resonant 
frequency. 
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Frequency Response from HP Spectrum Analyzer, 35 to 47.5 Hz 
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Figure 30 Damping of the Inertial Actuator/Platform System 


Frequency Response from HP Spectrum Analyzer 



Figure 31 Damping Calculation 
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E. CALIBRATION 

In order to use the LJC testbed, a calibration of the sensors and mirrors is 
necessary. Testing of the mirror’s motion about both axis and the sensor output is 
required to verify operability of the system. Each day, before testing, a calibration 
program is run to check the system. If at any time during testing, unusual results are 
obtained, the calibration program is run again to check system response. The program 
may be run from any directory, simply by typing the command “calibrate”. Once the 
calibration run is complete, the user may opt to keep the calculated gains for the system, 
or delete them if desired. 

To calibrate the system, a zero input data run is made. This shows how close the 
alignment is to having the beams “zeroed” on the detectors. Since the incident angle of 
the beam on the RFSM and detectors is critical to the ability of the system to control the 
beam, this zero input data run allows the user to align the folding mirrors to place the 
beam as close to zero as possible, at least within 0.1 mm, with the mirrors commanded at 
zero angle. The height of the X axis of the detectors, and the center of the RFSM, has 
been precisely placed at the height of the laser beam. By aligning the folding mirrors in 
this manner, as much of the nonlinearities (due to the angle of the beam) as possible is 
taken out of the system. Once that is complete, each mirror is commanded to rotate about 
its X and Y axis from the negative to the positive limits. With the laser beam on, this 
sweeps the beam across the PSDs in the horizontal (X axis) and vertical (Y axis) 
direction. Graphs are constructed using the detector output, allowing the user to check 
the response of both the detector and mirror. Step inputs are given to each axis of the 
mirrors, to check the step response of the system. Next, the mirrors center the beam on 
each of the PSDs, and calculate the input necessary to “zero” the beam. Using the data 
generated during the run of the detector response, gain of the mirrors, and the zero 
position, a 1 mm box is drawn centered on each detector. The user can easily verify 
proper operation by observing the boxes drawn by the mirrors. The figures below show 
the result for OT2, and compare the drawn boxes for a properly aligned system, vs. a 
poorly aligned system. The figure on the left is the three PSD outputs for a zero control 
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signal to the RFSM and DFSM prior to “drawing” the boxes. The difference in the OT2 
box and the OT3 box is due to the difference in distance between the target and feedback 
detector. The calibration program and typical output are provided in Appendix C. Figure 
33 also shows the “noise” in the beam from the thermal variations in the room and from 
the laser source itself 


zero input ncise, 00,012 and 013 r^wpcrt calibration box, OT2 and OT3 




time dstsnce, mm 

Figure 32 Poorly aligned Calibration Boxes 


zero input ncise, OT1 ,OT2 and OT3 


-u 


— 012 Y 
- OT3X 





Ntev'port calibration box, 0T2 and OT3 



dstance, mm 


Figure 33 Properly aligned Calibration Boxes 
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V. NUMERICAL SIMULATION OL THE CONTROL PROBLEM 


Adaptive control worked well for the situation in which the platform containing 
the control system was stationary and the disturbance was introduced only by the DFSM. 
Once we began to induce vibrations on the platform however, there were several 
problems that appeared. The flexibility of the platform and associated mirrors affected 
the reference signal, causing degradation in the adaptive system’s ability to remove the 
jitter in the beam. For example, the Gradient Adaptive Lattice (GAL) controller did not 
work well as frequency approached 50 Hz. For another, bias at the target was removed 
by the LMS controller, vice just passing this “DC signal” through. This bias, depending 
on the value, would severely impact the ability of the controller to remove even a single 
narrow band frequency. Although the mirrors and control system had been modeled in 
MATLAB, we did not see these effects in the simulation. In order to study these effects, 
we decided to devise a simple system, one that could be analytically solved and the 
simulations verified. The system should be straight forward, yet analogous to the control 
problem we faced. 

A. THREE MASS SPRING AND DAMPER SYSTEM 

A simple three mass spring and damper system was chosen, as shown below: 



Figure 34 Three Mass System 
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As can be seen from the diagram, the disturbance force goes “though” the control 
point, similar to the effect on the testbed, in that the “disturbed” light beam goes 
“through” the RFSM. For this simple system, the value of M, the mass of the one¬ 
dimensional blocks, are all the same. Likewise, K, the spring constant and C, the 
damping coefficient, are the same, simplifying the analytical calculations. The positions 
of each mass are X \, XI and X3 respectively. The control problem is to maintain mass 3 
at the zero position, regardless of the force exerted at mass 1, by using a control force at 
mass 2. Since we are dealing with discrete systems in the experiment on the testbed, we 
decided to choose the value of M and K so that the first fundamental frequency would be 
similar to that we see on the platform. In a discrete system, the control algorithm must be 
calculated within the sample time of the system. We therefore use the same sample rate 
for the three mass system as for the experiment, namely 2 kHz. 

B. ANALYTICAL SOLUTION 

Since the values for the spring constants are the same, and the values for the 
damping coefficients are the same, the system exhibits proportional damping, allowing us 
to use modal decomposition to solve the coupled differential equations of motion. The 
modal solution can be shown to be: 

q t (t) = 4 sin {cot + ft ) + e~^‘ (C u cos (a)j) + C 2i sin (coj)) 
i = 1,2,3 

where: 

t is time 

q t is the i th modal coordinate 
ft* is the i th natural frequency 
4 is the i ,h modal damping coefficient 
ft* is the forcing frequency, and 

ft> /( = , the i"' damped natural frequency 

Additionally, 
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where f. is the i tb modal forcing function amplitude. Setting the initial conditions to 
zero, the constants C\ and C 2 may be found: 

C u =~A i sin($) 

—A 1 (<W cos ($) + £.( 0 . sin(^.)) (5.4) 


Finally, mapping the modal coordinates to the coordinates of the system with the 
matrix of ortho-normal modes used to decouple the differential equations, we obtain: 


x(t) 


X\ 

X2 

X3 


_P]q(t) 


(5.5) 


where: 


<j>. are the vectors of the ortho-normal modes 


This system of equations can be used to solve for the motion of each of the 
masses, given an input force at mass 1 or mass 2 or a combination of both. The equations 
were programmed into MATLAB, to verify the accuracy of the MATLAB simulation of 
the three mass model. 
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C. NUMERICAL SIMULATION OL THE THREE MASS MODEL 

For use in numerical simulations in MATLAB, the three mass model is posed in a 
State Space system of equations. Using standard dynamics equations the model is: 
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(5.6) 


where Z is the disturbance input at mass 1 and U is the control input at mass 2. In order 

lb -Sec - 

to obtain natural frequencies near the platform frequencies, we set M = 20 —-- and 

in 

K = 2 x 10 7 —. These values give natural frequencies of 70.8, 198.5 and 286.8 Hz. C 
in 

may be set to obtain the damping desired. The State-Space system of equation (5.6) is 
discretized using a sample rate of 2 kHz. A comparison of this model’s output at mass 3 
to the analytical solution for the motion of mass 3 for a damping of 5% of critical for the 
lb see 

first mode (C = 4420—--) is shown in Figure 35, for a sinusoidal forcing function of 

in 

50 Hz. In the 3-mass model, the forcing function is implemented by the sine wave source 
block provided in the Digital Signal Processing Toolbox of MATLAB (DSP Sine Wave). 
The DSP Sine Wave is configured for multiple frequencies in the discrete mode, at a 
sample rate of 2 kHz. 
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x io' 3 comparison of analytical vs. numerical solution at 50 Hz input to mass 1 



Figure 35 Comparison of Analytical Model to the Numerical Simulation for a 50 Hz 

forcing frequency 

The Bode plot of the model is given below for 0% and 5% damping. 


mass 1 to mass 3 



Figure 36 Bode Plot for three mass model, input at mass 1 
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mass 2 to mass 3 



Figure 37 Bode Plot for three mass model, input at mass 2 


1. LMS Control of the Three Mass Model 

A simple LMS controller is developed using MATLAB/Simulink standard blocks, 
patterned after the LMS filter provided in the Digital Signal Processing (DSP) Toolbox of 
Simulink. The LMS filter in the DSP Toolbox does not have the capability to use a 
“filtered-x” version of the controller. The block diagrams below show the modification. 



Figure 38 LMS Controller Block Diagram 
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Figure 39 Coefficient Update Block Diagram 


The buffers in Figure 38 have length equal to the order of the filter. The overlap 
for each buffer is set to 1 less than the length, thereby producing a vector of delayed 
inputs at each time step. The leakage factor in Figure 39 may be used if there is 
insufficient spectral excitation of the algorithm, such as the case of a noiseless sinusoid 
[Kuo, 1996, p. 36]. In the studies conducted in this research, the leakage factor was not 
used and this gain was set to one. The Normalization block was used if the normalized 
version of the LMS (nLMS) algorithm was desired. This block divides the input vector 
by the Euclidian-Norm-Squared of the vector. The value “mu” in Figure 39 is the 
adaptation rate for the LMS algorithm or the value for a in equation (2.27) if the nLMS 
algorithm is used. Throughout the rest of this study, the nFXLMS algorithm is used, and 
a is replaced using the symbol //. The number of stages or order of the filter was 
determined by running simulations with increasing number of stages to determine the 
optimum number for the LMS controller. For this problem, we detennined that 16 stages 
is optimum. The reference signal is implemented using a second DSP Sine Wave. The 
reference signal amplitude and phase may then be varied to study the effects on the 
controller. 

a. Effect of Phase Shift of the Reference Signal 

In order to observe the effect of a phase shift of the reference signal, a 
small adaptation rate was selected. As noted earlier, small adaptation rates result in a 
smaller value of the misadjustment, and a better steady state response. This can be seen 
in the figure below. In Figure 40, two 60 second runs were made, one with a large [U and 
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one with a small value of //. The disturbance was a 30 Hz signal with a varying bias 
(described in section c below). The plot shows the effect of misadjustment. 


3-mass control problem at 5 % damping. 30 Hz dist freq 



Figure 40 Effect of misadjustment on steady state solution to the three mass system 

The adaptation rates selected for the three mass model were on the order 
of 1 x 10" 4 for the value of fu in order to observe the effect of phase on the reference 
signal. Runs were conducted, varying the forcing frequency from 10 to 170 Hz. The 
amplitude of the disturbance was adjusted to provide approximately 350 microns of 
motion at mass 3 for each disturbance frequency. The phase of the reference signal was 
varied from -3.1 to +3.1 radians at each frequency, and the MSE over the last one second 
of a 10 second run was plotted. The effect of the phase shift is easily seen in the 
following figures. 
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Minimum MSE vs phase for various freq. 



Figure 41 MSE vs phase shift of reference signal for 10 to 40 FIz 



phase, rad 


Figure 42 MSE vs phase shift of reference signal for 50 to 100 Hz 
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Minimum MSE vs phase for various freq. 



Figure 43 MSE vs phase shift of reference signal for 110 to 170 Hz 

The jump in the 110 Hz MSE result in Figure 43 is due to a slight change 
in the disturbance amplitude that the program caused by the simulation being run in 2 
simulations, from -3.1 to 0 radians and then from 0 to +3.1 radians. These figures clearly 
show that if the optimum phase shift for the reference signal can be used, the 
convergence is drastically effected for small values of //. 

b. Adaptive Delay Filter 

To find this optimum, and use it to increase convergence, an Adaptive 
Delay Filter (ADF) was developed. The ADF uses variable delays to increase the 
correlation between the phase of the reference signal and the phase of the error, as 
discussed in Chapter II. The block diagram of the ADF is shown below. 
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Figure 44 Adaptive Delay Filter 

In Figure 44, the reference signal is input at input port 1 (In 1). This 
signal is then delayed, by the block labeled Variable Fractional Delay and output at 
output port 1 (refsig delay). The delay is calculated using equation (2.41). The value of 
ju A used in this diagram is bounded by equation (2.42) and is set to 0.01 for the 

simulation runs. The constant “dels” is the estimate of the number of samples delay 
through the secondary plant, and the value “ddr” is normally set to 1 for the calculation 
involving the symmetric difference (see equation (2.41) and accompanying discussion). 
Since there is the competing effect between the LMS controller and the ADF of both 
driving the error to zero, one must give the ADF a good estimate of a starting delay, in 
order for it to have an effect. Otherwise, the error may be driven to zero before the ADF 
can adapt. Reviewing the figures above with respect to phase shift, it is noted that a lead 
angle will give the optimum effect. Therefore, the initial estimate used for the ADF is a 
delay equal to 90% of one period for the frequency of interest. This gives the ADF a 
reasonable starting point, and allows the optimum response. As an example of the effect 
of the ADF, two simulations are run on the three mass system. The first run uses the 
LMS controller alone, and the second, the LMS and ADF controller. 
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Figure 45 Comparison of LMS control with LMS+ADF 

Figure 45 clearly shows that use of the ADF improves convergence speed 
for this type of control problem. Misadjustment has not been affected. 

c. Effect of Bias on the Control System 

In using adaptive control to remove unwanted signals, such as in an 
electronic signal, the bias in the signal is normally removed prior to filtering. In the 
control of a laser beam with “noise”, there is usually a bias associated with the signal due 
to the fact that the beam does not hit the center of the target. Generally, this bias is 
removed by means of a compensator, and the adaptive filter/controller removes the noise. 
We noted during testing on the LJC testbed, that a small amount of bias added to the 
reference signal contributed significantly to the convergence of the algorithm. This effect 
was discussed in Chapter II, and is examined here using the three mass model. In 
observing the bias on the LJC testbed, a certain amount of randomness is evident. The 
figure below illustrates the laser beams x coordinate on the target sensor over 1.5 seconds 
with no movement of control mirrors: 
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Figure 46 Bias of Laser Beam on Target 

This variance is due to the effect of thermal gradients in the air 
surrounding the optical bench and platform, and due to the variance in the laser source. 
In order to accurately examine the bias effect on the three mass model, a similar varying 
bias was constructed using random noise and an averaging block in MATLAB. 


Bias 



Figure 47 Variable Bias simulated in MATLAB 
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The Adaptive Bias Filter (ABF) discussed in Chapter II is used to add bias 
to the reference signal. The figures below illustrate control with and without the ABF, 
for small values of // on the order of 1 x I O' 4 . The ABF was initialized using an estimate 
of C’ in equation (2.46) of 70% of the value of the mean bias present. Figure 48 shows 
the effect of the added bias. The position of mass three is centered at a value different 
from zero, corresponding to the bias in the disturbance signal. In Figure 49 this bias is 
removed. The Mean Square Error plot in Figure 50 clearly shows the advantage of 
adding bias to the reference signal. 
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Time Plot - 30 Hz disturbance control using LMS 
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Figure 49 Time Plot - 30 FIz disturbance control using LMS+ABF 


3-mass control problem at 5 % damping. 30 Hz dist freq 



Figure 50 MSE Plot - 30 Hz disturbance with varying bias 
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Unlike the ADF, the ABF does have an effect even at high values of //. 
Without any bias in the reference signal, the LMS controller will not “zero” the signal, 
resulting in a constant MSE, much higher than an LMS controller with bias added to the 
reference. This can be seen in the MSE plot for two simulations, one without bias and 
one with the ABF adding bias. 


3-mass control problem at 5 % damping. 30 Hz dist freq 



sec 


Figure 51 Comparison of LMS and LMS+ABF at high values of adaptation 

d. Effect of Anti-resonance Point on the Control System 
In Figure 37 an anti-resonance point can be seen at about 225 Hz. Control 
at this frequency would be expected to be not as effective as at other frequencies due to 
the loss of effectiveness of the motion at mass 2 affecting the motion at mass 3. This is 
borne out in Figure 52. At the anti-resonance point at 0.1% and to some extent 1% 
damping, control is not as effective as at the frequencies on either side of this point. This 
effect diminishes as the damping is increased. The effect is similar to that which is seen 
when an actuator is placed near a node for a prevalent mode shape in a structure. As the 
platform and control system did not exhibit an anti-resonance point, we were unable to 
verify this experimentally. 
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Figure 52 Effect of Anti-resonance point on LMS control 

2. Gradient Adaptive Lattice (GAL) Control of the Three Mass Model 

Another form of adaptive control is that which uses an adaptive lattice as 
discussed in Chapter II. The algorithm developed by Haykin [Haykin, 2002, p. 542] was 
modified for use with MATLAB/Simulink. The algorithm could not be implemented 
using standard Simulink blocks and was therefore coded for use with an S function. The 
C code for the S function is provided in Appendix D. Again, the optimum number of 
stages was determined in the same manner as the LMS filter discussed in section V.C.l 
above. 16 stages were used in the following numerical simulations. The reference signal 
was also provided by a separate DSP Sine Wave source block, in order to adjust bias, 
phase and amplitude of the signal. 

The GAL is similar to LMS control in that it is affected by phase differences and 
bias. However, GAL is much more sensitive than the LMS controller in these areas. The 
GAL algorithm uses an adaptation coefficient in a manner similar to the normalized LMS 
algorithm in that it varies with time, but instead of being normalized by the power of the 
reference signal, the GAL adaptation rate is normalized by the “power” in the prediction 
error. Thus, similar to the LMS we can evaluate slow and fast adaptation rates. We 
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cannot, however, compare values of adaptation rates between controllers, due to the 
difference in the normalization signals. Unlike the LMS algorithm, the effect on the 
misadjustment is slight for this type control problem, as can be seen in the following plot: 


GAL 3-mass control problem at 5 % damping. 40 Hz dist freq 



Figure 53 Misadjustment in GAL controller 
The speed of convergence is rapid for higher values of fi as can be seen below: 


GAL 3-mass control problem at 5 % damping. 40 Hz dist freq 



Figure 54 Speed of convergence for the GAL 
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a. Effect of Phase Shift of Reference Signal 

With the speed of convergence so fast, and the misadjustment slight for 
large values of jU for this problem, correlating the reference and error signals does not 
contribute appreciably to the solution of the control problem. However, there is an effect, 
similar to the LMS controller, for low values of //. Therefore, for problems with large 
eigenvalues of the input correlation matrix R, a low value of jU may be necessary and the 
use of a phase shift would aid in convergence (see equation (2.28)). We note that the 
amount of phase shift to produce an effect on the GAL convergence is far less than that 
necessary for the LMS convergence. The GAL algorithm is much more sensitive to the 
shift than the LMS, and if the phase is shifted too far, the algorithm becomes unstable. 


3-mass control problem at 5 % damping. 30 Hz dist freq 



Figure 55 Effect of Phase shift of reference signal on GAL with low values of // 
b. Effect of Bias on the GAL Control System 

Again, similar to the effect of bias on the LMS system, the GAL is 
affected by bias as well. A varying bias is applied to the disturbance signal as in the 
LMS problem. Bias is then added to the reference signal. Two effects are noted. First, 
the addition of bias to the reference signal allows the algorithm to converge faster, as well 
as center the third mass’ position at zero. Second, without bias added to the reference 
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signal, any increase in the value of /u results in an unstable controller. With bias added, 
/Li may be increased by an order of magnitude, resulting in rapid convergence and 
centering mass 3 at zero. 


3-mass control problem at 5 % damping. 30 Hz dist freq 



Figure 56 Effect of bias on the GAL controller 


3. Comparison of GAL and LMS 

Our goal is to rapidly correct the bias and narrow band frequency present in the 
disturbance. Our method of comparison will be the Mean Square Error (MSE). Our 
measure of effectiveness will be the convergence time, and the misadjustment. The 
misadjustment will be compared by observing the MSE for each controller at steady state. 
In the following simulations, each run will have a varying bias, as in Figure 47. We will 
start with the 30 FIz disturbance previously discussed, but then evaluate the controllers at 
frequencies above the first fundamental, as well as multiple frequencies. Both the LMS 
and GAL controllers will use 16 th order filters. 

a. Rate of Convergence and Misadjustment for a Single Frequency 

In order to compare the two controllers, the adaptation rate will be 
increased until the point just prior to instability. This will give us the fastest convergence 
for each controller. We will then compare the convergence rate and misadjustment. 
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Steady state will be determined by stability of the tap gains. Due to the emphasis on 
convergence at high values of fi , the ADF will not be used for the LMS controller. In 
these comparisons, the GAL controller will have the optimum bias added to the reference 
signal, and the LMS controller will use the ABF. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

sec 


Figure 57 Comparison of MSE for GAL and LMS, 3-mass system at 30 Hz 

As can be seen in the figure above, the GAL controller has a slightly faster 
convergence. The misadjustment is very nearly the same, with the mean value of the 
MSE for the LMS controller over the last 10 seconds of a 20 second run being 0.5 
micron larger than the GAL controller. A comparison of the steady state response is 
shown in Figure 58. A comparison of the Power Spectral Density (PSD) plot is shown in 
Figure 59. The PSD is calculated over the last 5 seconds of the 20 second simulation. 
From the PSD, we see that the LMS controller removes more of the 30 Hz disturbance 
than the GAL, but excites the first mode of the system as can be seen by the energy added 
at 70 Hz. 
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3-mass control problem at 5 % damping. 30 Hz dist freq 



Figure 58 Comparison of Steady State for LMS, GAL, 30 FIz disturbance 


PSD Plot 30 Hz dist freq 



Figure 59 Comparison of controller’s PSD response for 30 FIz 


Now we will compare controllers for a single frequency above the first 
fundamental of 70 FIz. The same simulation with varying bias will be run again at 87 Hz. 
An interesting phenomenon occurs above 70 Hz (the first phase shift between input and 
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output of the 3-mass system) for the GAL controller. It becomes unstable. In order to 
stabilize the output of the controller, the output of the GAL filter must be inverted and the 
adaptation rate significantly decreased. If the GAL controller is used to remove an 87 Hz 
electronic signal, (i.e. the disturbance does not go through a transfer function), it works as 
expected and no inversion is required. It is only when there is a transfer function 
between the disturbance and the output that provides a phase shift that the instability 
occurs. As discussed above, the GAL controller is very sensitive to phase. 


3-mass control problem at 5 % damping. 87 Hz dist freq 



sec 


Figure 60 


Comparison of MSE for GAL and LMS , 3-mass system at 87 Hz 
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3-mass control problem at 5 % damping. 87 Hz dist freq 



Comparison of Steady State for LMS, GAL, 87 Hz disturbance 
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Figure 61 








































































PSD Plot 87 Hz dist freq 



Figure 62 Comparison of controller’s PSD response for 87 FIz 

As can be seen from the above figures, the LMS controller converges 
faster and has lower misadjustment than the GAL controller at these frequencies. The 
LMS controller removes the narrowband disturbance, but still adds energy to the first 
mode, as in the 30 Hz case. 

b. Rate of Convergence and Misadjustment for Multiple 
Frequencies 

Our first comparison will be for two frequencies below the first 
fundamental. Both narrowband components have the same amplitude. Figure 63 shows 
that the LMS controller actually converges faster than the GAL. Additionally, the LMS 
controller has a better steady state response than the GAL. In order to minimize the 
instabilities in the GAL controller, low values of fi must be used, which leads to slower 
convergence. 
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Figure 63 


Figure 64 


3-mass control problem at 5 % damping. 30, 53 Hz dist freq 



Comparison of MSE for GAL and LMS ,3-mass system at 30 and 53 Hz 


3-mass control problem at 5 % damping. 30, 53 Hz dist freq 



Comparison of Steady State for LMS, GAL, 30 and 53 Hz disturbance 
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PSD Plot 30 53 Hz dist freq 



Figure 65 Comparison of controller’s PSD response for 30 and 53 FIz 


Our second simulation of multiple frequencies involves one narrowband 
frequency on either side of the first fundamental. In this case, the GAL is unstable except 
at very low adaptation rates, or when an IGRS is used. 


3-mass control problem at 5 % damping; 50 87 Hz dist freq 



Figure 66 Comparison of MSE for GAL and LMS ,3-mass system at 50 and 87 Hz 
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3-mass control problem at 5 % damping; 50 87 Hz dist freq 



Figure 67 Comparison of Steady State for LMS, GAL, 50 and 87 Hz disturbance 


PSD Plot 50 87 Hz dist freq 



Figure 68 Comparison of controller’s PSD response for 50 and 87 Hz 

As can be seen in Figure 66 through Figure 68, the GAL controller is not 
effective in controlling multiple frequencies, when one frequency is above the first phase 
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shift. Therefore, in systems in which multiple modes are excited, if some of the natural 
frequencies are on either side of a phase shift in the frequency response function, the 
GAL will have trouble controlling the system, depending on the relative amplitude of the 
narrowband peaks. 

c. The Addition of Random Noise to the Forcing Function 

Since we are also interested in the effect of random noise in the signal for 
the LJC testbed, random noise was inputted into the forcing function for the 3-mass 
model. The control of 200 Hz band-limited white noise in combination with two 
narrowband frequencies, 50 and 87 Hz as before, was simulated. For this case, we 
remove the variable bias in order to just examine the effect of random noise in the forcing 
function. Additionally, we use an IGRS for the GAL controller to correct the phase 
difference that occurs for the 87 Hz frequency. 


PSD Plot 50 87 Hz dist freq with 200 Hz band-limited white noise 



frequency, Hz 


Figure 69 Comparison of controllers PSD response for multiple frequencies and 

white noise 


Once again, the GAL is not as effective in controlling the two frequencies 
as the nFXLMS controller. 
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D. SUMMARY 

Simulation of the control problem using the 3-mass model allowed us to evaluate 
the effect of bias and phase shift on the reference signal. The simulation showed that a 
phase shift of the reference signal could increase the rate of convergence for an adaptive 
filter, without increasing the misadjustment. The addition of bias to the reference signal 
could correct a biased error signal, even one that is slowly varying. The simulation also 
shows that the GAL filter is very sensitive to the difference in phase between the 
reference signal and the error, resulting in instabilities unless the phase difference is 
corrected by means of an IGRS. 
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VI. DISTURBANCE REJECTION EXPERIMENTS ON THE LJC 

TESTBED 


The disturbance rejection experiments on the testbed were performed in three 
phases. The first phase involved the control of the laser beam with the platform 
stationary. This allowed us to evaluate the perfonnance of the control system under ideal 
conditions, as well as verify the operation and setup of the testbed. The second phase 
involved shaking the platform with discrete frequencies, while maintaining the DFSM at 
a constant position. This allowed us to see the effect of shaking the control system, and 
observe the control of vibratory motion by the RFSM. The third phase consisted of a 
randomly disturbed laser beam by the DFSM, and vibration of the platfonn. This would 
be what we would consider the on-orbit or on-station model. Our measure of 
effectiveness will be once again the Mean Square Error (MSE), determined by a running 
average filter of window 500 of the squared error. We will use the value of one standard 
deviation of the laser beam’s position about its mean to evaluate the “tightness” of the 
beam - how far it has spread out. We also measure the mean value of the beam’s 
position to determine how close to the center of the target the controller is placing the 
beam. 

A. LQR CONTROLLER FOR THE TESTBED 

In order to compare the classical control systems with the adaptive, an LQR 
controller was constructed using equations (4.1) and (4.2) as the state space system of 
equations. The linear quadratic optimal gains are detennined by using the command 
“LQRD” in MATLAB with a sample rate of 2 kHz, which calculates the required gains 
for the discrete system. The weighting matrix ( Q )for the states of the cost function (see 
equation (2.1) ) is an identity matrix, except for the values corresponding to the states for 
the voltage at the detector which are given a weight of 1000. The input weighting matrix 
(R ) is given a value of 0.1. The Linear Quadratic Regulator requires the measurement 
of each state in the state space system. In the testbed, we measure the voltage at the 
detector and the position of the RFSM, but not the velocity of the mirror in the RFSM, 
which is one of the states of the system. A Kalman estimator was constructed to estimate 
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the velocity of the mirror. The plant model in this state space system does not have an 
integrator. For the LQR to center the beam on target, an integrator was added to the 
model after the technique described in Ogata [Ogata, 2002, pp. 847-850]. 

The linear quadratic optimal gains calculated in this manner were then adjusted to 
obtain the tightest controlled beam possible. The program used to calculate the state- 
space system and the resulting output are provided in Appendix E. 

B. STATIONARY PLATFORM EXPERIMENTS 


Figure 8 is reprinted here for ease of reference. 



In the stationary platfonn experiments, the vibration isolation platfonn and the 
Newport Table are floated. The DFSM is used to inject a disturbance (narrowband, 
broadband or a combination) in the laser beam from the source. The Position Sensing 
Detector (PSD) labeled OT1 is used as a reference signal or as a signal to synchronize the 
Internally Generated Reference Signal (IGRS). The PSD labeled OT2 is the target for the 
beam. The power supply for the shaker is secured. The output from the PSDs are labeled 
X and Y axis for each PSD. The same controllers used for the 3-mass model are used for 
the testbed. The number of stages or order of the filter was determined by running 
experiments with increasing number of stages to detennine the optimum number for each 
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of the controllers. For the LJC testbed, 24 stages was optimum for the nFXLMS 
controller and 18 stages for the GAL controller. The “filter” for the nFXLMS controller 
is the state-space system of equations (4.1) and (4.2). The bias added to the IGRS for the 
nFXLMS controller is the mean value of the bias present at the OT1 sensor, and though 
bias is in the reference signal, it is not optimized. 

1. Periodic Disturbances- Phase Effect 

In the same manner as for the 3-mass model, the IGRS was phase shifted to 
determine the effect on the convergence of the nFXLMS controller for low adaptation 
rates. The experiment was run for 5 seconds, with the controller starting at 1.6 seconds 
into the run. This allowed time to record the disturbance at the target prior to control, in 
order to make comparisons before and after the controller was running. The MSE was 
calculated for the last 1 second of the run and used for comparison. The figures below 
provide the results for the Y axis of the PSD. The X axis is similar. The adaptation rate 
for the nFXLMS algorithm was 0.15 for the X-axis of the FSM. 


MSE vs Phase for various freq, 24 stage nFXLMS 



Figure 71 Effect of Phase shift of the IGRS for the testbed, low frequency 
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Min MSE vs Phase for various freq, 24 stage nFXLMS 



Figure 72 Effect of Phase shift of the IGRS for the testbed, high frequency 

As can be clearly seen, the phase shift has an effect on the convergence. In fact, 
for 50 through 70 Hz, it has a dramatic impact, improving the response by a factor of 5 
over the non-phase shifted IGRS for the 60 Hz disturbance. 

The GAL was also tested for the stationary platform case with periodic 
disturbance. The GAL controller uses the reference signal provided by OT1. The 
following figure compares the MSE for the GAL controller with the MSE for the 
nFXLMS controller, reference signal from OT1, and with the MSE for the nFXLMS 
controller using an IGRS and the optimum phase shift determined from Ligure 71 and 
Ligure 72. The MSE was detennined from the last one second of a 4.4 second run. 
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stable platform, MSE vs Freq 



Figure 73 MSE vs. Frequency for nFXLMS, phase-shifted nFXLMS, and GAL 

controllers, Stationary Platform 

Figure 73 shows that the GAL controller develops a lower MSE at 4.4 seconds 
than the nFXLMS controller using the same reference signal, for frequencies below 80 
Hz. However, the phase-shifted nFXLMS develops even lower MSE than the GAL 
above 30 Hz. 

In order to compare convergence rates, a single 50 Hz disturbance was injected by 
the DFSM for one axis, and the time vs. MSE plotted. The IGRS was modified by the 
ADF shown in Figure 44 for one of the runs. The ADF was initialized with a phase shift 
equal to 90% of the period of the disturbance signal. Additionally, a plot of the effect of 
using an IGRS with no phase shift is included for comparison. The figure shows that the 
GAL does reach its final value faster than the LMS filters, as is to be expected for this 
single frequency, however the GAL steady state value is higher than the LMS. The 
nFXLMS controller with ADF provides a 40% improvement in convergence rate over the 
nFXLMS using the IGRS alone, and a 50% improvement over the nFXLMS using OT1 
as the reference signal. 
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MSE for 50 Hz disturbance, stable platform 



Figure 74 Comparison of controllers for stationary platform 

2. The Effect of Bias on the Reference Signal for the Testbed 

During our first attempts to use an IGRS with the nFXLMS algorithm, we noted 
that the controller was far worse (in some cases adding energy to the disturbance) using 
the IGRS than when OT1 was used as a reference. The OT1 reference signal had a bias 
in it, due to the off-target position of the laser beam. Prior to the addition of the IGRS, 
we noted that depending on the location of the beam on OT1, we would obtain differing 
results. When we added a bias to the IGRS, the controller “behaved”, and removed the 
periodic disturbance as well as centered the beam on target. This effect caused the 
analytical investigation discussed in Chapter II, and resulted in the ABF of Figure 6. The 
figures below illustrate the effect of incorrect bias on the nFXLMS controller. 
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Figure 75 


Figure 76 


24 stage LMS Controller: Improvement: X -22.0274 %, Y 16.8826 % Mean X : 548 nm Y : -965 nm OT2 



50 Hz Stable Platform. Mean value of OT1X = 687p Mean Value of OT1Y = -1736(x 
Std Dev of error: X Axis input - 45.0343(1, Output - 54.9542p; Y Axis input - 48.4538p, Output - 40.2735(1 



50 Hz disturbance, Stationary Platform, using bias signal from OT1 


24 stage LMS Controller: Improvement: X 81.6211 %, Y 88.5686 % Mean X : -90 nm Y : 53 nm OT2 



50 Hz Stable Platform. Bias value for X IGRS = 26p Bias value for Y IGRS = -40p 


Std Dev of error: X Axis input - 45.2265p, Output - 8.3121 p; Y Axis input - 48.5737p, Output - 5.5526p 



50 Hz disturbance, Stationary Platform, using bias signal from ABF 
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Unlike the simulation of the 3-mass model, the nFXLMS controller on the LJC 
testbed would center the beam on the target regardless of whether bias was present or at 
the correct value in the reference signal. In the simulation of the 3-mass model, the 
controller would remove the noise, but leave the bias in the output, if no bias was present 
in the reference signal (see Figure 48). The effect of the bias in the testbed experiments 
causes a slower convergence to a steady state value, but the bias is removed. The reason 
for this is as yet unknown. It may be that the bias signal on the testbed is varying enough 
to cause the LMS algorithm to try and remove it, but this has not been verified to date. 
a. Proper Bias vs the Use of a Compensator 

The nFXLMS controller with proper bias was compared with a regular 
nFXLMS controller, no bias, using a parallel LQR as a compensator to determine if the 
addition of bias to the reference signal was better than using a compensator. A 50 and 87 
FIz signal was injected by the DFSM and the two controllers used to remove the error in 
the target signal. The PSD and MSE of the experiments are provided below. 


Stable Platform 2 freq (50, 87 Hz) PSD Plot - Y axis 



frequency, Hz 

Figure 77 PSD comparison of compensator and biased reference signal 
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MSE , stable platform 2 freq, 50 and 87 Hz 
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Figure 78 MSE comparison of compensator and biased reference signal 

From these figures, we see that the addition of bias to the reference signal 
results in a similar decrease in the PSD as using a straight compensator, and that the time 
constant for the biased reference signal is less than that for the parallel 
adaptive/compensator controller. 

3. Random Disturbances 

The DFSM was used to inject a random disturbance to the incoming laser beam. 
This random noise was a 200 Hz, band-limited white noise input to the DFSM, causing a 
corresponding random disturbance to the beam. This allowed us the ability to observe 
how the different controllers could handle random disturbances. 
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0 stage LQR Controller: Improvement: X 71.218 %, Y 65.1705 % Mean X : 236 nm Y :-274 nm OT2 



Std Dev of error: X Axis input - 40.6767p, Output - 11.7076p; Y Axis input - 34.6195(1, Output - 12.0578p 



Figure 79 LQR control of 200 Hz band-limited white noise 


Figure 79 shows that the LQR controller does a nice job of removing the noise, 
with a 71% improvement in the standard deviation of the beam pattern in the X position, 
and a 65% improvement in the Y position. The beam is centered on the detector, within 
the sensitivity of 500 nm of the sensor. But we must note that, due to the tight beam 
control desired, the gains are such that the integration to the steady state value requires 
about 6 seconds (see Figure 82). 

The figures below illustrate the nFXLMS and GAL controller response to the 
random disturbance, using the PSD signal from OT1 as the reference signal. The LMS 
algorithm does not do as well as the LQR, with only a 13% and 43% improvement in the 
X and Y position respectively. However, unlike the LQR, the nFXLMS and GAL 
controller reach their steady state value within 1 second (see Figure 82). 
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24 stage LMS Controller: Improvement: X 43.4785 %, Y 42.1462 % Mean X : 54 nm Y : 29 nm OT2 



Std Dev of error: X Axis input - 42.0776p, Output - 23.7829p; Y Axis input - 35.554p, Output - 20.5693p 



Figure 80 nFXLMS control of 200 Hz band-limited white noise 


Likewise, the GAL controller was tested with random noise, and the results are 
shown below. 


95 


































18 stage GAL Controller: Improvement: X 65.8207 %,Y 67.0166 % Mean X :-20 nm Y:-112nm OT2 



Std Dev of error: X Axis input - 44.2264p, Output - 15.1162^; Y Axis input - 37.2837|i, Output - 12.2974p 



Figure 81 GAL control of 200 Hz band-limited white noise 


The GAL controller exhibits the best control of the adaptive controllers for the 
stationary platform, random noise case, and approaches the LQR controller in 
effectiveness, as can be seen in the MSE plot below. The value (3 in the graphs below is 
the leakage factor discussed in Chapter II. 
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Figure 82 MSE comparison for random noise, stationary platform 

4. Multiple Frequencies 

The control of multiple frequencies for the stationary platform was attempted. 
The ADF system was not used for the nFXLMS due to the inability of the system to 
discriminate and phase shift multiple frequencies. We compared the LQR with the 
nFXLMS and GAL controllers (both using the OT1 output as the reference signal). The 
results for two frequencies injected by the DFSM are shown below. The nFXLMS does a 
better job of removing both frequencies than either LQR or GAL. 
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Stable Platform 50, 87 Hz PSD Plot 



Figure 83 Multiple frequency control - GAL and nFXLMS 

C. VIBRATING PLATFORM EXPERIMENTS 

Our next series of experiments involved maintaining the DFSM at a zero angle, 
and shaking the platform with the shaker mounted vertically (Y axis direction for the 
PSD). Due to the support structure of the platform and the mounting location, the 
platform moves in both the Y and X directions. The reference signal was synchronized 
with the disturbance using the either the signal to the shaker or the output from the 
accelerometer. 

1. Periodic Disturbance - Phase Effect 

The experiments conducted for the stationary platform configuration were 
repeated with the platform vibration providing the disturbance. The IGRS was again 
phase-shifted, and the results are shown below. Bias was added to the reference signal, 
equivalent to the bias seen at the reference PSD, OT1. The data was calculated from the 
last 1 second of a 4.4 second run. Similar to the stationary platform case, phase shift has 
an impact on the convergence properties of the LMS filter for this case. 
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MSE vs Phase for various freq, 24 stage nFXLMS 



Figure 84 Effect of Phase-shift on MSE for nFXLMS controller, vibrating platform, 

low frequency 


Min MSE vs Phase for various freq, 24 stage nFXLMS 



Figure 85 Effect of Phase-shift on MSE for nFXLMS controller, vibrating platform, 

high frequency 
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As in the stationary platform case, the GAL controller was tested at various 
frequencies and compared against the nFXLMS in its various configurations - with OT1 
as the reference signal and with the IGRS phase shifted to the optimum value. 


vibrating platform, MSE vs Freq 



Figure 86 MSE vs. Frequency for nFXLMS, phase-shifted nFXLMS, and GAL 

controllers, Vibrating Platform 

It is easy to see the effect of phase shifts on the GAL controller from Figure 86. 
From Figure 29, (reprinted below as Figure 87) we see that there is a phase difference 
between OT1Y and OT2Y starting at about 40 FIz, and a 180 degree difference at 80 FIz. 
Figure 86 shows that the GAL controller has an increase in MSE at 40 Hz, and actually 
adds energy to the disturbance at 80 Hz. 
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Y axis Frequency Response Function using 400 avgs. Sample Rate 2 kHz. Freq res 0.033333 Hz. NyFreq 1 kHz 




Figure 87 Frequency Response Function, OT1Y and OT2Y 


As in the stationary platform case, convergence rates are compared for a single 
vibration frequency of 50 Hz. Again, the nFXLMS using OT1 as the reference, the 
nFXLMS with the IGRS, and the nFXLMS with the IGRS and ADF combination are 
compared against the GAL controller. The figure below shows once again that the GAL 
controller converges to its steady state value the fastest (approximately 0.4 seconds for 
the situation). The nFXLMS with IGRS and ADF reaches the steady state value of the 
GAL controller output in approximately 2 seconds, and is 53% of the steady state value 
of the GAL controller. 
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MSE for 50 Hz disturbance, vibrating platform 



Figure 88 Comparison of controllers for vibrating platform at 50 Hz 

2. Periodic Disturbances - Bias Effect 

The vibrating platform experiments exhibit the identical response as the stationary 
platform with respect to the addition of bias to the reference signal. The addition of the 
proper bias to the reference signal allows for much faster convergence, with the same or 
better steady state MSE. 

3. Random Disturbances 

200 Hz band-limited white noise was injected into the shaker. The time domain 
results for each controller are shown below. These graphs show that for random 
vibrations, the LQR controller again works best. 
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LQR Controller: Improvement: X 21.9437 %, Y 68.9286 % Mean X : 119 nm Y : -267 nm OT2 



Std Dev of error: X Axis input - 21.0908|i, Output - 16.4627(i; Y Axis input - 49.712pi, Output - 15.4462(1 



Figure 89 LQR control of random vibrations 


24 stage LMS Controller: Improvement: X -28.0756 %, Y 18.738 % Mean X : -137 nm Y : -725 nm OT2 



fi x = 0.15 (i y = 0.05 


Std Dev of error: X Axis input - 16.3953|i, Output - 20.9984p; Y Axis input - 39.8097pi, Output - 32.3502|i 



time,sec 

Figure 90 LMS control of random vibrations 
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18 stage GAL Controller: Improvement: X -0.80188 %, Y 2.6487 % Mean X : 8 nm Y : -584 nm OT2 



Std Dev of error: X Axis input - 17.5581 p, Output - 17.6989p; Y Axis input - 32.0525p, Output - 31.2035p 



Figure 91 GAL control of random vibrations 

4. Multiple Frequencies 

The shaker was excited with multiple frequencies to evaluate the ability of the 
controllers to handle this type of disturbance. Again, 50 and 87 Hz were the input 
frequencies. As can be seen in the figure below, the GAL controller added energy to the 
87 Hz vibration, as well as the 100 Hz vibration from the 2 nd folding mirror. Both the 
nFXLMS and GAL controllers added energy to the resonance shaker frequency of 38 Hz. 
Based on the results of these experiments and the simulations on the 3-mass model, we 
determined that an IGRS was necessary for the GAL controller to work in a system in 
which a phase shift occurs between the reference signal and the error. 
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Vibrating Platform at 50, 87 Hz: PSD Plot 



Figure 92 PSD graph of two frequency control, vibrating platform 

D. COMBINATION OF VIBRATING PLATFORM AND EXTERNAL 
DISTURBANCES 

The final test of the LJC testbed involves the vibration of the platfonn at multiple 

frequencies and the DFSM injecting noise into the laser beam. This experiment is very 

close to the on-orbit or on-station condition expected for a laser relay platform. The 

uplink laser beam is corrupted by random disturbances due to thermals in the 

transmission path, or other causes. The platform itself is disturbed by on-board 

equipment that causes vibration of the structure, the total disturbance resulting in 

unwanted jitter. As discussed earlier, this jitter will result in a spreading of the beam, 

lowering the intensity at the target. Both the LMS and GAL controllers were configured 

with an IGRS to mitigate the effects of phase shifts between reference and error signals. 

Additionally, the LMS controller used the ABF filter. The reference signal for the GAL 

controller was provided a constant bias (optimized for the best response). Due to the 

complexities inherent in the GAL algorithm, the ABF filter could not be used with the 

GAL controller without overloading the CPU. The ADF was not used for either 

algorithm, as it does not contribute significantly to convergence at these adaptation rates. 
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1. Parallel LQR and nFXLMS Controllers 

As shown in section VLB.3 and VI.C.3 above, the LMS controller does not 
remove random noise well, but the LQR control is effective with this type of disturbance. 
By combining the two controllers in parallel, the rapid removal of bias and control of 
narrowband frequencies by the adaptive algorithms could enhance the removal of noise 
accomplished by the LQR controller. However, when combining the adaptive controller 
with the LQR or GAL controller, instabilities result when the LQR controller has an 
integrator. By setting the integration gain to zero, this instability is removed, and the 
adaptive controller then acts as a very fast integrator. 

2. The Vibration and Noise Experiments 

These experiments consist of two narrowband frequencies injected into either the 
DFSM or the shaker, one on either side of the phase shift between OT1 and OT2 for the 
vibrating platfonn case. 50 Hz and 87 Hz were chosen and used for both the DFSM 
narrowband disturbance and the shaker disturbance for comparison. In addition, the 
DFSM is used to inject 200 Hz band-limited white noise to simulate a random 
disturbance in the transmission path for the noise case. Each adaptive controller was 
used separately and in parallel with the LQR controller. The experiment was performed 
in four cases, covering the stationary and vibrating platform multiple frequency and noise 
situations. In each case, the Power Spectral Density plot and Mean Square Error plot are 
compared for each of the control schemes. The Power Spectral Densities are calculated 
using the last 4 seconds of data of a 10 second run, by Welch’s method. The data is 
windowed using a 2048 length Hamming Window, and the Fast Fourier Transform (FFT) 
performed using a 1024 length FFT with a 50% overlap. Results are tabulated for each 
case. In the following tables, the Input Jitter is the recorded sensor readings at the target 
in the first second of the run, prior to the controller starting. The standard deviation of 
each axis for the input is recorded. The Controlled Beam is the one second average of the 
standard deviation of the beams position for that axis, at the end of the 10 second run. 
The Mean value of the beam position is the average of the beam’s X or Y coordinates 
over the last one second of the 10 second run, in nano-meters (nm). Note that the 
minimum sensitivity of the PSD is ±500 nm. The amount that each of the two tones is 
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reduced is listed (+ indicates energy added to that frequency). The Total MSE is the 
combined MSE for each axis, using the averaging filter discussed previously. The Time 
Constant is the time from controller cut-on until the Input MSE is reduced by a factor of 
1/e. The PSD plots for each of the cases are for the Y axis only, since the major 
disturbance axis for the vibration case is in the Y direction. The X axis results are 
provided in Appendix F. 

a. Case 1: Stationary Platform Narrowband Disturbance from 
DFSM 

For this case, the platform remains stationary. 50 and 87 Hz narrowband 
frequencies are summed and injected into the DFSM to provide the disturbance. The 
PSD, MSE and tabulated results follow. 


Stable Platform 2 freq PSD Plot - Y axis 



Figure 93 PSD Case 1, DFSM narrowband disturbance 
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Power Spectral Density, dB/Hz 


Stable Platform 2 freq PSD Plot - Y axis Combo controllers 



Figure 94 PSD Case 1, Combination Controllers 


Total MSE Stable Platform 



Figure 95 MSE Case 1 
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Total MSE , Stable Platform Combo controllers 



Figure 96 MSE Case 1, Combination Controllers 


Controller 

LQR nFXLMS GAL 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

45.4 

49.8 

45.8 

49.9 

48.1 

52.1 

Controlled Beam, Std Dev., 
microns 

10.9 

16.0 

9.9 

6.6 

21.1 

16.2 

No. of stages/order 

n/a 

n/a 

24 

24 

18 

18 

Adaptation rate, fj. 

n/a 

n/a 

0.10 

0.05 

0.10 

0.06 

% reduction in jitter 

76.0 

68.0 

78.3 

86.9 

56.1 

68.9 

Mean value of beam 
position, nm 

60 

-595 

-695 

-38 

165 

-61 

db reduction in PSD of 50 
Hz 

-20.4 

-17.0 

-30.6 

-30.5 

-14.4 

-27.3 

db reduction in PSD of 87 
Hz 

-14.1 

-12.0 

-16.6 

-20.8 

-0.1 

-2.0 

Total MSE at 10 seconds, 
micron 2 

372.0 

142.4 

680.1 

Time constant, secs 

0.71 

0.27 

0.26 


Table 3 Case 1 results 
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Controller 

i LQR+nFXLMS LQR+GAL ] 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

45.3 

49.9 

47.9 

52.0 

Controlled Beam, Std Dev., 
microns 

10.1 

14.1 

19.7 

15.5 

No. of stages/order 

24 

24 

18 

18 

Adaptation rate, fj. 

0.10 

0.05 

0.10 

0.06 

% reduction in jitter 

77.6 

71.8 

58.8 

70.1 

Mean value of beam 
position, nm 

54 

52 

67 

-165 

db reduction in PSD of 50 
Hz 

-38.4 

-39.2 

-24.5 

-21.3 

db reduction in PSD of 87 
Hz 

-24 

-23.3 

-11.4 

-9.1 

Total MSE at 10 seconds, 
micron 2 

299.0 

596.0 

Time constant, secs 

0.16 

0.16 


Table 4 Case 1, Combination Controller results 


Case 1 shows that the nFXLMS controller works best in removing the 
stationary platform DFSM narrow band disturbance. The addition of the LQR controller 
adds to the time required for the controller to approach its steady state value, and results 
in a higher MSE than for the controller acting alone. This higher MSE for the 
combination controller is due to the addition of energy by the LQR controller at the 
higher frequencies, around 600 Hz, as can be seen in the PSD for this case (Figure 97) 
when expanded to include frequencies to 800 Hz. The LQR controller alone exhibits a 
similar high frequency response. However, the combination controller does remove the 
narrowband component better than the nFXLMS controller alone, and centers the beam 
better on the target. It is noted here, that in each of the 4 cases, the steady state MSE 
value for the GAL or GAL/LQR controller is rapidly reached, independent of the type of 
disturbance. For this particular geometry, the GAL controller reaches its steady state 
value, in each of the cases, in less than 0.4 seconds. As discussed in section II.C, the 
GAL filter is expected to rapidly converge to its steady state value, and this is borne out 
by these experiments. 
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Stable Platform 2 freq PSD Plot - Y axis 
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Figure 97 Case 1 Comparison of nFXLMS and nFXLMS/LQR 

b. Case 2: Stationary Platform Narrowband Disturbance Plus 
Noise from the DFSM 

In this experiment, the platfonn is again stationary. 50 and 87 Hz 
narrowband frequencies are summed and injected with 200 Hz band limited white noise 
into the DFSM to provide the disturbance. The PSD, MSE and tabulated results follow. 
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Figure 98 
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PSD Case 2, DFSM narrowband disturbance plus 200 Hz noise 

Stable Platform + noise 2 freq PSD Plot - Y axis Combo controllers 
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Figure 99 PSD Case 2, Combination Controllers 
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Total MSE Stable Platform + noise 



Figure 100 MSE Case 2 


Total MSE , Stable Platform + noise Combo controllers 



Figure 101 MSE Case 2, Combination Controllers 
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Controller 

LQR nFXLMS GAL 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

62.3 

63.7 

63.1 

63.2 

66.6 

65.2 

Controlled Beam, Std Dev., 
microns 

12.3 

16.3 

57.0 

52.3 

25.3 

20.4 

No. of stages/order 

n/a 

n/a 

24 

24 

18 

18 

Adaptation rate, // 

n/a 

n/a 

0.1 

0.05 

0.1 

0.06 

% reduction in jitter 

80.3 

74.5 

9.7 

17.2 

62 

68.7 

Mean value of beam 
position, nm 

87 

-515 

-1915 

3321 

-13 

-319 

db reduction in PSD of 50 
Hz 

-20.7 

-18 

-9.9 

-12.9 

-17.8 

-17.7 

db reduction in PSD of 87 
Hz 

-14.3 

-12.3 

-0.8 

-1.1 

-0.7 

-3.7 

Total MSE at 10 seconds, 
micron 2 

416.0 

5979.0 

1165.0 

Time constant, secs 

0.74 

0.3 

0.25 S 


Table 5 Case 2 Results 


Controller 

i LQR+nFXLMS LQR+GAL | 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

63.2 

63.4 

65.7 

67.6 

Controlled Beam, Std Dev., 
microns 

12.3 

15.3 

22.8 

17.5 

No. of stages/order 

24 

24 

18 

18 

Adaptation rate, fl 

0.1 

0.05 

0.1 

0.06 

% reduction in jitter 

80.6 

75.8 

65.3 

74.2 

Mean value of beam 
position, nm 

-63 

-23 

20 

-167 

db reduction in PSD of 50 
Hz 

-34.8 

-37.8 

-25 

-21.5 

db reduction in PSD of 87 
Hz 

-20 

-21.1 

-12.6 

-10.5 

Total MSE at 10 seconds, 
micron 2 

386.0 

917.0 

Time constant, secs 

0.16 

0.16 


Table 6 Case 2 Results, Combination Controllers 


In this case, we see the detrimental effect of random noise on the 
nFXLMS controller. The random noise is almost entirely passed through the controller. 
This can be expected since the reference signal has none of the random component. Even 
if the reference signal has the random component, as discussed in section VI.B.3 above, 
the nFXLMS controller does not respond well. Surprisingly, the GAL controller 
outperforms the nFXLMS controller, again with only the narrowband components in its 
reference signal. Although there is an increase in energy for this controller above 150 
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Hz, the random component is almost entirely rejected in the region below 100 Hz, a 
result not expected based on the analysis of the GAL controller in section II.C. The 
combination nFXLMS/LQR controller works best for this case. The narrowband 
components are removed by the LMS algorithm and the random component is removed 
by the LQR controller. Although the LQR adds energy above 500 Hz as discussed 
above, the removal of the random component in the range below 150 Hz more than 
makes up for this in the MSE. The parallel nFXLMS/LQR controller responds much 
faster than the LQR alone, as can be seen by a time constant of 0.16 versus 0.74 seconds. 

c. Case 3: Narrowband Disturbance from the Vibrating Platform 
In the Case 3 experiments, we use the inertial actuator to vibrate the 
platform at the same two narrowband frequencies used in Case 1 and 2 above and 
compare the results. In this case, the amplitude of the narrowband frequency in the X 
axis is less, due to the mounting configuration of the inertial actuator in the Y direction. 
The adaptation rate for the X-axis nFXLMS controller (X-axis of the mirror) was lowered 
to 0.05 from 0.1 in order to improve performance for the vibrating platform case. 


Vibrating Platform 2 freq PSD Plot - Y axis 



Figure 102 PSD for Case 3: Vibrating Platform Disturbance 
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Vibrating Platform 2 freq PSD Plot - Y axis Combo controllers 



Figure 103 PSD for Case 3, Combination Controllers 


Total MSE Vibrating Platform 



Figure 104 MSE for Case 3 
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Total MSE , Vibrating Platform Combo controllers 



Figure 105 MSE for Case 3, Combination Controllers 


Controller 

LQR nFXLMS GAL 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

26.4 

53.4 

24.5 

52.2 

24.3 

53.0 

Controlled Beam, Std Dev., 
microns 

9.9 

15.5 

12.0 

7.8 

19.4 

19.3 

No. of stages/order 

n/a 

n/a 

24 

24 

18 

18 I 

Adaptation rate, [X 

n/a 

n/a 

0.05 

0.05 

0.1 

0.06 

% reduction in jitter 

62.4 

71.1 

51.0 

85.1 

20.5 

63.6 

Mean value of beam 
position, nm 

18 

364 

-70 

318 

9 

-26 

db reduction in PSD of 50 
Hz 

-14.6 

-14.3 

-24.1 

-31.2 

-15.9 

-28.4 

db reduction in PSD of 87 
Hz 

-14.8 

-11.0 

-2.2 

-21.6 

-4.2 

-5.3 

Total MSE at 10 seconds, 
micron 2 

337.8 

203.7 

848.7 

Time constant, secs 

0.72 

0.26 

0.25 


Table 7 Case 3 results 
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Controller 

LQR+nFXLMS LQR+GAL | 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

27.3 

57.6 

24.1 

53.2 

Controlled Beam, Std Dev., 
microns 

9.8 

11.4 

19.2 

17 

No. of stages/order 

24 

24 

18 

18 

Adaptation rate, /! 



0.1 

0.06 

% reduction in jitter 

64.1 

80.1 

20.6 

68.1 

Mean value of beam 
position, nm 

20 

318 

1 

-125 

db reduction in PSD of 50 
Hz 

-29.5 

-37.7 

-19.5 

-21 

db reduction in PSD of 87 
Hz 

-16.0 

-28.0 

-13.3 

-9.8 

Total MSE at 10 seconds, 
micron 2 

225.0 

681.1 

Time constant, secs 

0.34 

0.16 


Table 8 Case 3 Results, Combination Controllers 


The addition of the IGRS to the GAL controller has mitigated the effect of 
the phase shift between reference signal and error as can be seen from the decrease in the 
energy at 87 Hz. As discussed in section VI.C.4, the GAL controller added energy for 
frequencies above the phase shift between OT1 and OT2. Here, the IGRS effectively 
shifts the phase 180 degrees for the 87 Hz reference signal, allowing the GAL filter to 
adapt to lower the energy in this frequency. For Case 3, the nFXLMS controller 
performs nearly as well as the combination nFXLMS and LQR controller. The 
combination controller is better in removing the narrowband components than the 
nFXLMS controller alone as in Case 3. 

d. Case 4: Narrowband Disturbance from the Vibrating Platform 
Plus 200 Hz Noise from the DFSM 

In the final case, 200 Hz band limited white noise from the DFSM is 
added to the disturbance from the inertial actuator. In order to graphically illustrate the 
case and its results, some time domain plots are provided for the best controllers. First, to 
demonstrate the effect of the disturbances on the laser beam, a two dimensional view of 
the experiment beam pattern is shown below. The graph on the right shows the effect of 
the 50 and 87 Hz vibration of the platform on the beam. The graph on the left shows the 
effect of the 50 and 87 Hz vibration, plus the random disturbance from the DFSM. Each 
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asterisk represents the beam’s position on the detector at 2 msec intervals. The blue 
circle in the diagram is drawn at the geometric center of the pattern of the laser beam, and 
the radius of the circle is one standard deviation of the pattern over 0.1 seconds. Each 
graph is a snapshot of the beam’s position over a 0.1 second interval. 
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Figure 106 Two Dimensional Beam Pattern 
The PSD and MSE results are shown and tabulated below. 


Vibrating Platform + noise 2 freq PSD Plot - Y axis 
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Figure 107 PSD for Case 4, Vibrating Platform plus 200 Hz noise from the DFSM 
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Vibrating Platform + noise 2 freq PSD Plot - Y axis Combo controllers 



Figure 108 PSD for Case 4, Combination Controllers 


Total MSE Vibrating Platform + noise 



Figure 109 MSE for Case 4 
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Total MSE , Vibrating Platform + noise Combo controllers 



Figure 110 MSE for Case 4, Combination Controllers 


Controller 

LQR nFXLMS GAL 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

47.7 

68.2 

47.2 

63.9 

51.0 

66.7 

Controlled Beam, Std Dev., 
microns 

12.4 

18.2 

43.3 

46.4 

25.5 

25.8 

No. of stages/order 

n/a 

n/a 

24 

24 

18 

18 

Adaptation rate, /z 

n/a 

n/a 

0.05 

0.05 

0.1 

0.06 

% reduction in jitter 

73.9 

73.3 

8.3 

27.3 

50.0 

61.4 

Mean value of beam 
position, nm 

165 

-763 

342 

-3063 

-84 

-373 

db reduction in PSD of 50 
Hz 

-14.1 

-15.1 

-5.0 

-16.1 

-7.5 

-11.9 

db reduction in PSD of 87 
Hz 

-13.3 

-11.2 

+1.6 

-7.5 

-0.5 

-4.3 

Total MSE at 10 seconds, 
micron 2 

484.8 

4054.9 

1444.7 

Time constant, secs 

0.73 

0.26 

0.25 


Table 9 Case 4 Results 
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Controller 

| LQR+nFXLMS LQR+GAL J 

Control Mirror Axis 

X-axis 

Y-axis 

X-axis 

Y-axis 

Input Jitter, Std. Dev, 
microns 

44.9 

65.5 

49.2 

66.4 

Controlled Beam, Std Dev., 
microns 

12.5 

15.3 

23.2 

19.2 

No. of stages/order 

24 

24 

18 

18 

Adaptation rate, /J, 


0.05 

0.1 

0.06 

% reduction in jitter 

72.2 

76.7 

52.8 


Mean value of beam 
position, nm 

2 

57 

-24 

-136 

db reduction in PSD of 50 
Hz 

-23.7 

-38.7 

-18.5 

-21.6 

db reduction in PSD of 87 
Hz 

-13.1 

-10.4 

-13.4 

-9.8 

Total MSE at 10 seconds, 
micron 2 

389.2 

936.5 

Time constant, secs 

0.18 

0.16 


Table 10 Case 4 Results, Combination Controllers 


The time domain data for the nFXLMS/LQR combination controller is 
shown below. These results indicate a 72% to 76% reduction in the standard deviation of 
the input jitter, and the mean value of the beam pattern is well within the measurement 
sensitivity of the PSD (± 500 mn). The GAL/LQR controller showed a slightly less 
improvement, 52% to 71% also with the center of the beam pattern within the 
measurement sensitivity. 
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24 stage MQG Controller: Improvement: X 72.2362 %, Y 76.6715 % Mean X : 2 nm Y : 57 nm OT2 



Std Dev of error: X Axis input - 44.9558p, Output -12.4815p; Y Axis input - 65.4481 p, Output - 15.268p 



time,sec 

Figure 111 Time Domain results (nFXLMS/LQR) for Case 4 


The 2-D views of the nFXLMS/LQR results are provided below. The 
target is the inner black circle, with a radius of 25 microns. The outer circle has a radius 
of 50 microns. 

Band-Limited White Noise @200Hz plus 50, 87 Hz vibration @l50p- 2 msec asterisk, 0.1 sec snapshot 



Figure 112 2-D view of the nFXLMS/LQR controller effect on optical beam 
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The total MSE vs. time graphs are most revealing for this case. In Figure 
110 we see the effect of using the adaptive nFXLMS controller in parallel with the FQR 
controller. By adding the proper bias to the reference signal using the ABF, the 
nFXFMS/FQR controller is able to rapidly zero the signal, and remove the noise. In 
combining the two controllers, we get the speed of the adaptive controller in centering the 
beam and removing the disturbance narrowband frequencies, and the ability of the FQR 
controller to remove noise, together resulting in a “tight” beam centered on the target. 


124 



VII CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE 

RESEARCH 


A. CONLUSIONS 

As part of this research to develop improved techniques for the control of Optical 
Beam Jitter, a Laser Jitter Control Testbed was designed and built at the Spacecraft 
Research and Design Center, Naval Postgraduate School to validate the new techniques 
experimentally. This Testbed is unique in its ability to experimentally test algorithms for 
the control of jitter in an optical beam. The testbed has the capability to vibrate the Fast 
Steering Mirror and associated optics used to correct the disturbances to the beam in a 
controlled manner, allowing the investigation of different algorithms in suppressing jitter. 

This research has shown that the phase difference between the reference signal 
and the error signal is important in correcting a narrowband disturbance. This phase 
difference develops since a physical structure is being excited by an inertial actuator or 
some other vibratory motion. As a result of the excitation, phase differences occur 
between the reference signal and error signal since these signals are measured at different 
locations on the structure. By correlating the reference and error signal, improved 
convergence rates may be obtained, without increasing the adaptation rate. Since the 
adaptation rate can be lower, we have increased the convergence rate without increasing 
misadjustment. 

The Gradient Adaptive Lattice filter is particularly susceptible to this phase 
difference. For the GAL filter to control vibrations in a structure, the disturbance must be 
below the point at which this phase difference occurs, or the reference signal must be 
phase-shifted to compensate. Thus the selection and modification of the reference signal 
and an understanding of the physical structure and its frequency response is necessary to 
properly optimize the GAL filter. The GAL filter does show promise in its ability to 
correct band-limited white noise that is not structure related - for example, random 
disturbances that occur as the beam transits the atmosphere. 
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We have also found that the addition of the correct amount of bias to the reference 
signal can have the affect of a compensator, allowing the adaptive filter to not only 
correct the disturbance, but rapidly remove any bias present in the error. This 
compensation method is faster than applying an integrator to the error signal in an 
adaptive system in which an LQR compensator is used. The LQR compensator can then 
be used without an integrator, specifically tuned to remove random noise. This 
compensator can be placed in parallel with the LMS controller using bias in the reference 
signal to control a “colored noise” situation, in which there is bias in the error. The 
resulting combination controller is faster and more effective at removing jitter than either 
the adaptive system or the LQR controller alone. 

In addition to the general conclusions above, several improvements in optical 
beam jitter control techniques have been made during the course of this research. In 
particular: 

• An Adaptive Delay Filter was developed and implemented to correlate the error 
and reference signal for the LMS controller. This filter, designed originally for 
the measurement of time delays in a sampled data system, is used to adaptively 
delay the reference signal to better correlate the signal with the error. 

• The use of an internally generated reference signal, phase shifted to compensate 
for the difference between the reference measurement point and the error, was 
developed and implemented to correct the phase difference problem described 
above for the Gradient Adaptive Lattice filter 

• An Adaptive Bias Filter has been developed and implemented to adaptively add 
bias to the reference signal for adaptive filters. The ABF allows rapid correction 
of bias errors by an adaptive filter without the use of a compensator. 

• A parallel LMS/LQR controller has been developed and implemented to remove 
bias in the error, correct narrowband disturbances from a vibrating platfonn 
supporting the control system, and control noise in an optical beam. 

B. RECOMMENDATIONS FOR FUTURE RESEARCH 

As a result of the work perfonned during the course of this research, several areas 
have arisen that require further study. These include: 

• The GAL filter derived for use in the controller for the testbed does not have a 
provision for the use of a “Filtered-X” version similar to the LMS algorithm. Kuo 
(1996) has developed a Filtered-X version of a Lattice filter. This version should 
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be adapted for use in the testbed, and may result in better performance due to the 
addition of the filtered reference signal. 

• The Gradient Based Lattice Filter developed by Chen and Gibson [Chen, 2001], 
may prove to be superior the GAL filter. This filter was used by Boelitz, et al 
[Boelitz, 2003] in the ABLE ACE experiment discussed in Section II. 
Preliminary simulations run on the mathematical model of the LJC testbed using 
this type filter show that it may outperform both the GAL and the nFXLMS filter 
in its ability to handle both random and narrowband disturbances. Modifying this 
filter for use on the testbed and subsequent testing will confirm its usefulness for 
controlling the laser targeting problem. Although not used on the testbed, the C- 
code is provided in Appendix E for use in MATLAB/Simulink. 

• It is still not fully understood, why the bias present in the error signal is removed 
when no bias is provided in the reference signal. Simulations and analysis 
indicate that the bias should remain, and the noise and/or narrowband 
disturbances should be removed. In the LJC testbed, the presence of bias in the 
error signal drastically affects the ability of the controller to remove the unwanted 
disturbance, yet the bias itself is removed. Perhaps it is the variability of this bias, 
but that must still be proved. 
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APPENDIX A: FAST STEERING MIRROR DATA 


A. NEWPORT FAST STEERING MIRROR 


Number of Axes 

2 (tip-tilt) 

Angular Range from ± 10 V 

±26.2 mrad (±1.5°), Mechanical'" 

Resolution 

<! grad mis. Mechanical 111 

Repeatability 

<3 grad mis, Mechanical 111 

Accuracy From ±26,2 mrad, 25°G ,U) 

<0.262 mrad (0.015°). Mechanical" 1 

Linearity From ±26.2 mrad, 25°C^ 1,2 * 

<1.0% 

Peak Angular Velocity 

>2.5 radVsec, Mechanical* 11 

Peak Angular Acceleration 

>900 Tad/sec*, Mechanical" 1 

Closed-Loop Amplitude Bandwidth* 7 ’(-3 dB) 

i 

>550 H/ at 100 (trad amplitude 

Optional: >750 Hz at I0U grad amplitude 

Closed-Loop Phase Bandwidth 121 (60° lag) 

>300 Hz at 100 prad amplitude 

Gain Margin 

>10 dB 

Phase Margin 

245° 

Response Flatness 1 ’ 1 

Peaking <3 dB 

Small Angle Step Response for Steps < 250 grad. 
Mechanical 11 ' 


Rise Time to 90% 

<3.5 msec 

Settling Time to 5% 

<5.0 msec 

Large Angle Step Response for 26.2 mrad Steps, 
Mechanical* 1 


Rise Time to 90% 

<3.5 msec 

Settling Time to 5% 

<12 msec 

Cross-axis Coupling, Static 

<0.1% 

Cross-axis Coupling, Dynamic 

<0.1% 

Powered Null Offset (Closed-Loop) 

±1 mrad. Mechanical* 1 ' 

Noise Equivalent Angle (1 Hz to 10 kHz) 

<3 grad rms 

Resolution of 3x>cal Position Sensor 

<0.5 pirad 

Quiescent Power at FSM Assembly 

<5 W at any angle ±26.2 mrad 

Operating Temperature Range* 2 ' 

Otu 50°C (32 to 122°F) 

Storage Temperature Range 

-20 to 55%’ (-4 to 131 °F) 

Warm-up Time for Mirror Stability 121 at 25°C 

<10 minutes 

Mirror Thermal Drift 121 

<5 prad/ c C, Mechanical*" 

Optical Axis Location without base 

1.5 in. (38.1 mm) High, Centered Lefr-to- 
Right 

Mass with base 

1,1 lb (0.5 kg) 

Envelope without base 

3.0 x 3.0 * 2.3 in. [w x h. x dl 
(76.2 x 76.2 < 58.7 mm) 


9.S ft (3 in) 


Table 11 Newport Fast Steering Mirror Data 
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Typical FSM-CD100 Controller/Driver Specifications 


Command Input and DIT Output 

Analog, ±10 V =±26.2 mrad 

Peak Operating Power to Mirror 

30 W 

Continuous Maximum Operating Power to 
Mirror 

15 W 

Thermal Protection 

60°C at mirror coil 

Current Protection 

3 A 

Operating Temperature Range 

0 to 50°C (32 to I22°F) 

Storage Temperature Range 

-20 to 55°C (-4 to 131°F) 

Power i 

100-240 V, 50/60 Hz, 2.5 A 

Mass 

10.7 lb (4.9 kg) 

Envelope (without rack mount flanges) 

17 x 3.5 x 11.59 in. [w * h x d] 

(431.8 x 88.9 x 294.4 mm) 


Table 12 Newport Fast Steering Mirror Controller/Driver Specifications 



Figure 3: Safe Operating Area for FSM-200 Series fast steering mirrors 


Figure 113 Newport Fast Steering Mirror Safe Operating Area 
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B. BAKER ADAPTIVE OPTICS FAST STEERING MIRROR 


LIGHT FORCE ONE KEY SPECIFICATIONS 


SMALL SIGNAL RESPONSE 


LARGE SIGNAL RESPONSE 


cm ‘ SfiRWw:'ic' 'SijumvfcM iftms' 


(■tit "# J i> : u Ttrryv 


ifcfcVB MJUW 


IANDWIDTH VS, 


RISETIME 


VS. STEP amplitude 


NOISE EQUIVALENT ANGLE; 
MIRROR: 


<20 nanoRadlanslll 

25.4mm dia., 6.35mm thick. <1/10A @633nm 


MECHANICAL 

POWER 


6" Square x 2 " thick w/ample 1/4-20 mounting holes , wt.=5.5lbs 
25 Watts max. 


Contact Info: Jeffrey T. Baker, BAKER ADAPTIVE OPTICS 

1155 Sandia Dr., Bosque Farms, NM 87068, JTBaker@flash.net, 505-259-6628 phone, 505-869-2725 fa 

Figure 114 Baker Adaptive Optics - “Light Force One” 
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ONTRAC POSITION SENSING DETECTOR 

Model 2L10SP 


Active Responsitivity Dark Current Noise Capacitance 

Area @940nm nA Current PF@15V 

(mm) TYP TYP MAX TYP MAX 

MAX 


2L10SP 

10.0x10.0 

0.63 

100 500 

1.3 

2.5 

90 110 

Model 

Rise 

Reverse 

Detector 

Thermal 

Position 


Time 

10-90% 

15V 

TYP MAX 

Bias V 

MIN TYP 
MAX 

Resistance 

(K-ohm) 

MIN TYP MAX 

Drift 

ppm/C* 

TYP MAX 

Non- 

Linearity 

+/-% 
TYP MAX 

2L10SP 

.40 .80 

5 15 20 

7 10 16 

40 

200 

0.3 0.8 


Table 13 OnTrac Position Sensing Detector Specifications 



Figure 115 


OnTrac Position Sensing Module Diagram 



























































D. 


CSA INERTIAL ACTUATOR 



Engineering, Inc. 


Actuator number: SA5 - SN 116 

Test date: 2 December 2003 
Current drive amplifier 

Mechanical interface: single Vi-20, both ends 
Gender mate, and mating cable included 

Resistance. trnomJeijjoeca.ture.)vj!b.Utarls.V„.' , .5_Qhms_ 

Inductance (room temperature with leads): 1.03 mH 

Maximum continuous drive current, 0-peak: 4.0A 



CSA Engineering Inc., 2565 Leghorn St., Mountain View, CA, 94043-1616, USA 
www.csaengincering.com 650 210 9000 


Figure 116 CSA Inertial Actuator Response Characteristics 
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APPENDIX B: SOFTWARE VERSIONS 


Software 

Version 

Function 

MATLAB 

6.5 R13 

Interface Control Computer 
with xPC Target Computer 

SIMULINK 

5.0 R13 

Model and controller 
Software interface 

RealTime Workshop 

5.0 R13 

Interface Control Computer 
with xPC Target Computer 

xPC Target 

2.0R13+ 

Interface Experiment with 
Control Computer 

dSPACE 

Release 3.3 

Interface Experiment with 
Disturbance Computer 

Microsoft Visual C++ 

6.0 

Compiler 

DSP Blockset 

5.0 R13.0.1 

MATLAB Toolbox 

Control System Toolbox 

5.2 R13.0.1 

MATLAB Toolbox 

Windows XP 

5.1 Build 2600 Service 

Pack 1 

Control and Disturbance 
Computer Operating 

System 


Table 14 Software Versions 
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APPENDIX C: CALIBRATION PROGRAM 


A. CALIBRATION PROGRAM 

% file to calibrate mirrors and detectors 

% Mirror gains will be placed in a file in the working directory called 
% "mgains.mat" 

% Output figures are in directory "C:\Watkins\MATLAB Files\output figs\ 
% This Program uses model calibration_control 


Sample Time for model: do not set less than 0.0005 


Ts=0.0005; 


Initiate input/output variables 


otlx=l; 

otly=2; 

ot2x=5; ot2y=6; ot3x=7; 

ot3y=8; 

nxin=10; 

nyin=9; 

bxin=4; byin=3; 


signl=-l; 

const=l; 

ramp=l; 


tstng=-l; 

step=2; 

step_time : 

=0.25; 


% 





% distances in mm from mirrors to detectors 


dl=1585; % Dist from DFSM Mirror to OT1 

d2=370; % Dist from RFSM Mirror to OT3 

d3=1245; % Dist from RFSM Mirror to OT2 
dll=1370; % Dist from DFSM Mirror to OT3 


Initiate switches 


-s 


x_t e s t=ramp; 

y_test=ramp; 

bx=0; 

by=0; 

b_test_x=const; 

b_test_y=const; 

b_x_zero=0; 

b_y_zero=0; 

n_x_zero=0; 

n_y_zero=0; 

n_test_x=const; 

n_test_y=const; 

rampstep=step; 

finval=l; 

bak_x=l; 

bak_y=l; 

newxl=l; 

newyl=l; 

newx2=l; 

newy2=l; 

ramp_slope=10; 

ramp_initial=-10; 




tol=0.001; %tolerance for determining zero position 

% do not set less than 0.001 


Initiate Program 


disp ('energize Laser, OT301 amps, mirrors. Float Table if desired'); 
disp('allow 20 min stabilization time for accurate readings'); 
curr_dir=cd; 

cd('C:\Watkins\Matlab Files\Calibration programs'); 

disp('If Calibration model is open, it must be manually compiled before running') 
disp('If model is open and not built, exit this program by typing x and enter') 
replyl = input('otherwise press Enter to continue ','s'); 

if ~isempty(replyl) 

disp('Manually Build calibration_control model, or exit model, and restart 
calibrate'); 

error('exiting program'); 

end 


Check connections 


% is connection with target working? 
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if ~strcmp(xpctargetping, 'success') 

error('Connection with target cannot be established'); 
end 

% is calibration_control already open? 
systems = find_system('type', 'block_diagram'); 
if isempty(strmatch('calibration_control', systems, 'exact')) 
mdlOpen = 0; 

load(xpc,'calibration_control'); 
else 

mdlOpen = 1; 
end 


Obtain values for parameters from model 


tg=xpc; 

b_test_x_par=getparamid(tg,'Switching/b test x','Value'); 
b_test_y_par=getparamid(tg,'Switching/b test y','Value'); 

rampstep_par=getparamid(tg,'Test Signals/rampstep','Value'); 
ramp_slope_par=getparamid(tg,'Test Signals/Ramp/Step','After'); 
ramp_initial_par=getparamid(tg,'Test Signals/Ramp/Constant1','Value'); 


finval_par=getparamid(tg,'Test Signals/Stepl',' 

After' ) ; 

n_test_x_par=getparamid(tg, 
n_test_y_par=getparamid(tg, 

'Switching/n 
'Switching/n 

test x' 
test y' 

,'Value'); 
,'Value'); 

b_x_zero_par=getparamid(tg, 
b_y_zero_par=getparamid(tg, 

'Switching/b 
'Switching/b 

x zero' 
y zero' 

,'Value'); 
,'Value'); 

n_x_zero_par=getparamid(tg, 
n_y_zero_par=getparamid(tg, 

'Switching/n 
'Switching/n 

x zero' 
y zero' 

,'Value'); 
,'Value'); 


x_test_par=getparamid(tg,'Test Signals/zero','Value'); 
y_test_par=getparamid(tg,'Test Signals/zerol','Value'); 

bak_x_par=getparamid(tg,'bak_x','Gain'); 
bak_y_par=getparamid(tg, 'bak_y', 'Gain'); 

newx_par=getparamid(tg,'new_x','Gain'); 
newy_par=getparamid(tg,'new_y','Gain'); 

clc 


Set switches in model for zero input values 
and start model to obtain and plot data 


set_switches 

+tg 

pause (1.5); 

-tg 

tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'zero input', 'NumberTitle', 'off'); 

plot(tt,oo(:,otlx),tt,oo(:,otly),tt,oo(:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3 
y)),grid,zoom 

xlabel('time') 
ylabel('amplitude, V') 

legend COT1 X','OTl Y','OT2 X','OT2 Y','OT3 X','OT3 Y') 
title ('zero input noise, OTl,OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\zero_input_noise.fig') 
avgx=mean(oo(1000:2000, otlx)) ; 
avgy=mean(oo(1000:2000, otly)) ; 

disp('Zero value should be + or - 0.1V. Is the zero Satisfactory?'); 

reply = input('Enter n to manually adjust, otherwise, hit Enter to continue: ','s'); 


Loop for adjustment of zero if necessary 
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if ~isempty(reply) 

tg.StopTime=9999; 

+tg 

xpcscope 

disp('set mirrors to center positions on 0T1, 0T2 and 0T3,then Enter') 
disp('OTl is scope 1, 0T2 is scope 2, 0T3 is scope 3') 
disp('Load Scope "align_scopes" if desired') 

user_entry = input('hit Enter when aligned and ready to proceed'); 

-tg 

b_test_x=const; b_test_y=const; 

rampstep=step; 

set_switches 

+tg 

pause (1.5); 

-tg 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'zero input, adjusted', 'NumberTitle', 'off'); 

plot (tt,oo (:,otlx),tt,oo(:,otly),tt,oo(:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3y)) 
,grid,zoom 

axis([-inf inf -0.1 0.1]) 

xlabel('time') 

ylabel('amplitude, V') 

legend COTl X','0T1 Y','OT2 X','OT2 Y','OT3 X','OT3 Y') 
title ('zero input noise, OTl,OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\zero_input_noise.fig') 
avgx=mean(oo(1000:2000,otlx)); 
avgy=mean(oo(1000:2000, otly)) ; 

end 


Begin test of mirrors and obtain calibration data 


% Set the values for a ramp test of DFSM 
rampstep=ramp; 

ramp_slope=10; ramp_initial=-10; 

b_test_x=tstng; b_test_y=const; 


% ramp DFSM in theta x to get OT limits 


set_switches 

+tg 

pause(2.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'Baker Ramp X, OT1', 'NumberTitle', 'off'); 

plot(tt,oo(:,otlx),tt,oo(:,otly)),grid,zoom 

xlabel('time') 

ylabel('amplitude, V') 

legend('OT1 X','0T1 Y') 

title('Baker ramp response,theta X, OT1') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Baker ramp theta X, OTl.fig') 


Calculate min and max values from PSD 


nlymin=min(oo(250:4000,otly)); 
nlymax=max(oo(250:4 000, otly)) ; 
flagt=0; 

for n=250:round(1/Ts); 

tahd=n+round(0.05/Ts) ; 

if (~flagt)&((oo(tahd,otly)-oo(n,otly))>0.2) ; 
Bvlx_min=oo(tahd,bxin) ; 

Vdetly_min=oo(tahd, otly) ; 
break; 


139 









end; 


end 

for n=round(1.5/Ts):4000; 
tahd=n+round(0.05/Ts); 
if (oo(tahd,otly)-oo(n,otly)<0.2) ; 
Bvlx_max=oo(n,bxin); 
Vdetly_max=oo(n,otly); 
break; 

end 

end 


% ramp DFSM in theta y to get OT limits 


b_test_x=const; b_test_y=tstng; 

set_switches 

+tg 

pause(2.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'Baker Ramp Y, OT1', 'NumberTitle', 'off'); 

plot (tt,oo(:,otlx),tt,oo(:,otly)),grid,zoom 

xlabel('time') 

ylabel('amplitude, V') 

legend('OT1 X','0T1 Y') 

title('Baker ramp response,theta Y, OT1') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Baker ramp theta Y, OTl.fig' 


Calculate min and max values from PSD 


nlxmin=min(oo(250:4000,otlx)); 
nlxmax=max(oo(250:4000, otlx)) ; 
for n=250:round(1/Ts) ; 

tahd=n+round(0.05/Ts) ; 

if (~flagt)&((oo(tahd,otlx)-oo(n, otlx))>0.2) ; 
Bvly_min=oo(tahd, byin) ; 

Vdetlx_min=oo(tahd, otlx) ; 
break; 

end; 

end 

for n=round(1.5/Ts):4000; 
tahd=n+round(0.05/Ts); 
if (oo(tahd,otlx)-oo(n,otlx)<0.2) ; 

Bvly_max=oo(n,byin); 

Vdetlx_max=oo(n, otlx) ; 
break; 

end 

end 


% ramp RFSM in theta x to get OT limits 


% set switches for RFSM test 
b_test_x=const; b_test_y=const; 

n_test_x=tstng; n_test_y=const; 

set_switches 

+tg 

pause(2.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 
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set(figure, 'Name', 'Newport Ramp response', 'NumberTitle', 'off'); 

plot(tt,oo(:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3y)),grid,zoom 
xlabel('time') 
ylabel('amplitude, V') 

legend C0T2 X','0T2 Y','0T3 X','0T3 Y') 

title('Newport ramp response,theta X, 0T2 and 0T3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport ramp response,theta X, 0T2 and 
0T3.fig') 


Calculate min and max values from PSD 


n3ymax=max(oo(250:4000,ot3y)); 
n3ymin=min(oo(250:4000,ot3y)); 
n2ymax=max(oo(250:4000, ot2y)) ; 
n2ymin=min(oo(250:4000,ot2y)); 
for n=250:round(1/Ts) ; 

tahd=n+round(0.05/Ts) ; 

if (~flagt)&((oo(tahd,ot3y)-oo(n, ot3y))>0.2) ; 
Nv3x_min=oo(tahd, nxin) ; 

Vdet3y_min=oo(tahd,ot3y); 
break; 


end; 


end 

for n=round(1.5/Ts):4000; 
tahd=n+round(0.05/Ts); 
if (oo(tahd,ot3y)-oo(n,ot3y)<0.2) ; 
Nv3x_max=oo(n, nxin) ; 
Vdet3y_max=oo(n, ot3y) ; 
break; 


end 


end 


% ramp RFSM in theta y to get OT limits 


n_test_x=const; n_test_y=tstng; 

set_switches 

+tg 

pause(2.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'Newport Ramp Y, OT2,OT3', 'NumberTitle', 'off'); 
plot (tt,oo (:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3y)),grid,zoom 
xlabel('time') 
ylabel('amplitude, V') 

legend('OT2 X','OT2 Y','OT3 X','OT3 Y') 

title('Newport ramp response,theta Y, OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport ramp response,theta Y, OT2 and 
OT3.fig') 


Calculate min and max values from PSD 


n3xmax=max(oo(250:4000, ot3x)) ; 
n3xmin=min(oo(250:4000, ot3x) ) ; 
n2xmax=max(oo(250:4000,ot2x)); 
n2xmin=min(oo(250:4000, ot2x)) ; 
for n=250:round(1/Ts) ; 

tahd=n+round(0.05/Ts) ; 

if (~flagt)&((oo(tahd,ot3x)-oo(n,ot3x))>0.2) ; 
Nv3y_min=oo(tahd,nyin) ; 

Vdet3x_min=oo(tahd, ot3x) ; 
break; 

end; 

end 

for n=round(1.5/Ts):4000; 
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tahd=n+round(0.05/Ts) ; 
if (oo(tahd,ot3x)-oo(n,ot3x)<0.2) ; 
Nv3y_max=oo(n, nyin) ; 
Vdet3x_max=oo(n,ot3x); 
break; 

end 

end 


% calibration calcs mm/V on OT1/OT2 


otlxcal=10/(nlxmax-nlxmin); otlycal=10/(nlymax-nlymin); %mm/V 

ot2xcal=10/(n2xmax-n2xmin); ot2ycal=10/(n2ymax-n2ymin); 

ot3xcal=10/(n3xmax-n3xmin); ot3ycal=10/(n3ymax-n3ymin); 




NGmy=(Vdet3x_max-Vdet3x_min)*ot3xcal/d2/2*1000/(Nv3y_max-Nv3y_min); 
Y axis 

% 

mrad/V, 

Newport 

NGmx=(Vdet3y_max-Vdet3y_min)*ot3ycal/d2/2*1000/(Nv3x_max-Nv3x_min); 
X axis 

% 

mrad/V, 

Newport 

BGmx=(Vdetly_max-Vdetly_min)*otlycal/dl/2*1000/(Bvlx_max-Bvlx_min); 
axis 

% 

mrad/V, 

Baker X 

BGmy=(Vdetlx_max-Vdetlx_min)*otlxcal/dl/2*1000/(Bvly_max-Bvly_min); 
axis 

% 

mrad/V, 

Baker Y 


otcal=l; 


% calibration calcs mm/V on OT1/OT2 if ramp does not work 


% calc is voltage from max limits on OT. 

% for PSM 10-2 with OT 301 it is +/- 10 V or 20V total 
% across detector. Use if edge voltage of detector in doubt 
% clear flag otcal if desired to use these calcs 

if -otcal; 
calc=20; 

otlxcal=10/(calc);otlycal=10/(calc) ; 
ot2xcal=10/(calc);ot2ycal=10/(calc) ; 
ot3xcal=10/(calc);ot3ycal=10/(calc) ; 

end 


% step DFSM mirror in theta x 


rampstep=step;b_test_x=tstng; 

set_switches 

+tg 

pause(1.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 


% determine gain, V/mm 
% assume small angle approximation 


initl=mean(oo(10:4 00, otly)) ; 
finall=mean(oo(2000:2500,otly)); 
stepsizeyl=final1-init1; 
bak_xl=finval/(stepsizeyl*otlycal) ; 
bak_x=bak_xl; 

init2=mean(oo(10:400, ot2y)) ; 
final2=mean(oo(2000:2500, ot2y)) ; 
stepsizey2=final2-init2; 
bak_x2=finval/(stepsizey2*ot2ycal) ; 

init3=mean(oo(10:4 00, ot3y)) ; 
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final3=mean(oo(2000:2500,ot3y)); 

stepsizey3=final3-init3; 

bak_x3=finval/(stepsizey3*ot3ycal); 

set(figure, 'Name', 'Baker step X', 'NumberTitle', 'off'); 

plot(tt,oo(:,otlx),tt,oo(:,otly),tt,oo(:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3 

y)),grid,zoom 

xlabel('time') 

ylabel('amplitude, V') 

legend COT1 X','OTl Y','OT2 X','OT2 Y','OT3 X','OT3 Y') 
title('Baker step response,theta X, OTl,OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Baker step response,theta X, OTl,OT2 
and OT3.fig') 


% step DFSM mirror in theta y 


b_test_x=const; b_test_y=tstng; 

finval=l; 

set_switches 

+tg 

pause(1.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

initl=mean(oo(10:400, otlx) ) ; 
finall=mean(oo(2000:2500, otlx)) ; 
stepsizexl=final1-init1; 
bak_yl=finval/(stepsizexl*otlxcal) ; 
bak_y=bak_y 1; 

init2=mean(oo(10:400, ot2x)) ; 
final2=mean(oo(2000:2500, ot2x)) ; 
stepsizex2=final2-init2; 
bak_y2=finval/(stepsizex2*ot2xcal); 

init3=mean(oo(10:4 00, ot3x)) ; 
final3=mean(oo(2000:2500,ot3x)); 
stepsizex3=final3-init3; 
bak_y3=finval/(stepsizex3*ot3xcal) ; 

set(figure, 'Name', 'Baker Step Y', 'NumberTitle', 'off'); 

plot(tt, oo(:, otlx) , tt, oo(:,otly),tt,oo(:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3 

y)),grid,zoom 

xlabel('time') 

ylabel('amplitude, V') 

legend COT1 X','OTl Y','OT2 X','OT2 Y','OT3 X','OT3 Y') 
title('Baker step response,theta Y, OTl,OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Baker step response,theta X, OTl,OT2 
and OT3.fig') 


% find DFSM zero on OT1 


valx_sig=getsignalid(tg, 'OT output Scopes/On Trac/MM-32 /pi'); 

valy_sig=getsignalid(tg, 'OT output Scopes/On Trac/MM-32 /p2'); 

b_test_x=const; b_test_y=const; 

set_switches 

clear valx valy 

ycnt=0;xcnt=0; 

+tg 

for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
while abs(avgx)>tol 
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b_y_zero=b_y_zero-(avgx/abs(stepsizexl)); 
set (tg,b_y_zero_par,b_y_zero) ; 
for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
xcnt=xcnt+l; 
if xcnt>10 
break 

end 

end 

for i=l:5 

valy(i)=get(tg,valy_sig); 

end 

avgy= m ean(valy); 
while abs(avgy)>tol 

b_x_zero=b_x_zero-(avgy/abs(stepsizeyl)); 
set (tg,b_x_zero_par,b_x_zero); 
for i=l:5 

valy(i)=get(tg,valy_sig) ; 

end 

avgy= m ean(valy); 
ycnt=ycnt+l; 
if ycnt>10 
break 

end 

end 


% check b_y_zero again 


for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
xcnt=0; 

while abs(avgx)>tol 

b_y_zero=b_y_zero-(avgx/abs(stepsizexl)); 
set(tg,b_y_zero_par,b_y_zero); 
for i=l:5 

valx(i)=get(tg,valx_sig) ; 

end 

avgx=mean(valx); 
xcnt=xcnt+l; 
if xcnt>10 
break 

end 

end 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'Baker zero 0T1', 'NumberTitle', 'off'); 

plot(tt,oo(:,otlx),tt,oo(:,otly)),grid,zoom 

axis([-inf,inf,-1,1]) 

xlabel('time') 

ylabel('amplitude, V') 

legend('0T1 X','0T1 Y') 

title ('Baker zero, 0T1') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Baker zero.fig') 
b_y _ z e r o_o 11 =b_y_ zero; b_x_ z e r o_o 11 =b_x_ zero; 


% step RFSM mirror in theta x 


n_t e s t_x=t stng; 
rampstep=step; 
x_test=ramp; 


n_test_y=const ; 
finval=.5; 
y_test=ramp; 
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set_switches 

+tg 

pause(1.5); 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

init3=mean(oo(10:4 00, ot3y)) ; 
final3=mean(oo(1000:2500, ot3y)) ; 
init2=mean(oo(10:4 00, ot2y)) ; 
final2=mean(oo(1000:2500, ot2y)) ; 
stepsizey3=final3-init3; 
stepsizey2=final2-init2; 

newx3c=finval/(stepsizey3*ot3ycal);newx2c=finval/(stepsizey2*ot2ycal); 

set(figure, 'Name', 'Newport Step X', 'NumberTitle', 'off'); 

plot(tt,oo(:,ot2x),tt,oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3y)),grid,zoom 

xlabel('time') 

ylabel('amplitude, V') 

legend('OT2 X','OT2 Y','OT3 X','OT3 Y') 

title('Newport 0.5V step response,theta X, OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport 0.5V step response,theta X 
OT2 and OT3.fig') 


% step RFSM mirror in theta y 


n_test_x=const; n_test_y=tstng; 
finval=.5; 

set_switches 

+tg 

pause(1.5); 

-tg 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

init3=mean(oo(10:4 00, ot3x) ) ; 
final3=mean(oo(1000:2500, ot3x)) ; 
init2=mean(oo(10:400, ot2x)) ; 
final2=mean(oo(1000:2500,ot2x)); 
stepsizex3=final3-init3; 
stepsizex2=final2-init2; 

newy3c=finval/(stepsizex3*ot3xcal);newy2c=finval/(stepsizex2*ot2xcal); 

set(figure, 'Name', 'Newport Step Y', 'NumberTitle', 'off'); 

plot(tt, oo(:, ot2x) , tt, oo(:,ot2y),tt,oo(:,ot3x),tt,oo(:,ot3y)),grid,zoom 

xlabel('time') 

ylabel('amplitude, V') 

legend('OT2 X','OT2 Y','OT3 X','OT3 Y') 

title('Newport 0.5V step response,theta Y, OT2 and OT3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport 0.5V step response,theta Y 
OT2 and OT3.fig') 


% find RFSM zero on OT3 


b_test_y=const; b_test_x=const; 
n_test_y=const; n_test_x=const; 

fct=0.1; % Value for sensitivity of zero algorithm steps 

set_switches 

valx_sig=getsignalid(tg, 'OT output Scopes/On Trac/MM-32 /p6'); 
valy_sig=getsignalid(tg, 'OT output Scopes/On Trac/MM-32 /p5'); 
clear valx valy 
ycnt=0;xcnt=0; 
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+tg 

for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
while abs(avgx)>tol 

n_y_zero=n_y_zero-fct*(avgx/abs(stepsizex3)*finval); 
set (tg,n_y_zero_par,n_y_zero); 
for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
xcnt=xcnt+l; 
if xcnt>10 
break 

end 

end 

for i=l:5 

valy(i)=get(tg,valy_sig); 

end 

avgy= m ean(valy); 
while abs(avgy)>tol 

n_x_zero=n_x_zero-fct*(avgy/abs(stepsizey3)*finval); 
set (tg,n_x_zero_par,n_x_zero); 
for i=l:5 

valy(i)=get(tg,valy_sig); 

end 

avgy= m ean(valy); 
ycnt=ycnt+l; 
if ycnt>10 
break 

end 

end 


%check n_y_zero again 


for i=l:5 

valx(i)=get(tg,valx_sig) ; 

end 

avgx=mean(valx); 
xcnt=0; 

while abs(avgx)>tol 

n_y_zero=n_y_zero-fct*(avgx/abs(stepsizex3)*finval) ; 
set (tg, n_y_zero_par, n_y_zero) ; 
for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
xcnt=xcnt+l; 
if xcnt>10 
break 

end 

end 

-tg 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'Newport Zero 0T3', 'NumberTitle', 'off'); 

plot (tt,oo(:,ot3x),tt,oo(:,ot3y)),grid,zoom 

axis([-inf,inf,-1,1]) 

xlabel('time') 

ylabel('amplitude, V') 

legend('0T3 X','0T3 Y') 

title ('Newport zero, 0T3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport zero 0T3.fig') 


hold RFSM cal values for 0T3 


146 







n_y_ zero_ot 3=n_y_ zero; 
n_x_zero_ot 3=n_x_zero; 


% find RFSM zero on OT2 


finval=0.5; 
n_y_zero=0; 
b_test_y=const; 
n_test_y=const; 


n_x_zero=0; 
b_test_x=const; 
n_test_x=const; 


fctl=0.1; fct2=0.05; % Value for sensitivity of zero algorithm steps 

set_switches 

valx_sig=getsignalid(tg, 'OT output Scopes/On Trac/MM-32 /p3'); 
valy_sig=getsignalid(tg, 'OT output Scopes/On Trac/MM-32 /p4'); 
clear valx valy 
ycnt=0;xcnt=0; 

+tg 

for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
while abs(avgx)>tol 

n_y_zero=n_y_zero-fct1*(avgx*finval/stepsizex2); 
set (tg, n_y_zero_par, n_y_zero) ; 
for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
xcnt=xcnt+l; 
if xcnt>30 
break 

end 

end 

for i=l:5 

valy(i)=get(tg,valy_sig); 

end 

avgy=mean(valy); 
while abs(avgy)>tol 

n_x_zero=n_x_zero-fct1*(avgy*finval/stepsizey2); 
set (tg,n_x_zero_par,n_x_zero); 
for i=l:5 

valy(i)=get(tg,valy_sig); 

end 

avgy=mean(valy); 
ycnt=ycnt+l; 
if ycnt>30 
break 

end 

end 


%check n_y_zero again 


for i=l:5 

valx(i)=get(tg,valx_sig) ; 

end 

avgx=mean(valx); 
xcnt=0; 

while abs(avgx)>tol 

n_y_zero=n_y_zero-f ct2*(avgx*finval/stepsizex2) ; 
set(tg,n_y_zero_par, n_y_zero) ; 
for i=l:5 

valx(i)=get(tg,valx_sig); 

end 

avgx=mean(valx); 
xcnt=xcnt+l; 
if xcnt>30 
break 
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end 


end 

-tg 

clc 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

set(figure, 'Name', 'Newport Zero OT2', 'NumberTitle', 'off'); 

plot (tt,oo(:,ot2x),tt,oo(:,ot2y)),grid,zoom 

axis([-inf,inf,-1,1]) 

xlabel('time') 

ylabel('amplitude, V') 

legend('0T2 X','0T2 Y') 

title('Newport zero, 0T2') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport zero 0T2.fig') 

n_y_zero_ot2=n_y_zero; 

n_x_z e r o_ot 2 =n_x_z e r o; 


% draw boxes 

% set box size in mm: xdist or ydist is 1/2 of a side 


xdist=l; ydist=l; 


% draw box on 0T1 using DFSM 


bak_x=bak_xl*xdist; 
x_test=signl; 
b_t e s t_x=t s t ng; 
n_test_x=const; 
b_x_zero=b_x_zero_ot1; 
b_y_zero=b_y_zero_ot 1; 


bak_y=bak_yl*ydist; 
y_test=signl; 
b_test_y=tstng; 
n_test_y=const ; 


set_switches 

+tg 

pause (5.1); 
-tg 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 


% convert volts from sensor to mm 


oo(:,otlx)=oo(:,otlx)*otlxcal; oo(:,otly)=oo(:,otly)*otlycal; 

set(figure, 'Name', 'Baker Box 0T1', 'NumberTitle', 'off'); 

plot(oo(:,otlx),oo(:,otly)),grid,zoom 

axis([-5 5 -5 5]) 

xlabel('distance, mm') 

ylabel('distance, mm') 

legend('0T1') 

title ('Baker calibration box, 0T1') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Baker calibration box, OTl.fig' 


% draw box on 0T2 and 0T3 using RFSM 


newxl=newx2c*xdist; 
x_test=signl; 
b_test_x=const; 
n_test_x=tstng; 
n_x_zero=n_x_zero_ot2; 

set_switches 

+tg 

pause (5.1); 

-tg 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 


newyl=newy2c*ydist; 
y_test=signl; 
b_test_y=const; 
n_test_y=tstng; 
n_y_zero=n_y_zero_ot2; 
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% convert volts from sensor to mm 


oo(:,ot3x)=oo(:,ot3x)*ot3xcal; oo(:,ot3y)=oo(:,ot3y)*ot3ycal; 
oo(:,ot2x)=oo(:,ot2x)*ot2xcal; oo(:,ot2y)=oo(:,ot2y)*ot2ycal; 
set(figure, 'Name', 'Newport Box OT2,3', 'NumberTitle', 'off'); 

plot(oo(:,ot2x),oo(:,ot2y),oo(:,ot3x),oo(:,ot3y)),grid,zoom 
axis([-5 5 -5 5]) 
xlabel('distance, mm') 
ylabel('distance, mm') 
legend('0T2','0T3') 

title('Newport calibration box, 0T2 and 0T3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport calibration box, 0T2 and 

0T3.fig') 


apply offsets for 0T3 and draw 0T3 box 


n_y_zero=n_y_zero_ot 3; n_x_zero=n_x_zero_ot 3; 

newxl=newx3c*xdist; newyl=newy3c*ydist; 

set_switches 

+tg 

pause (5); 

-tg 

clear tt oo 
tt=tg.Time; 
oo=tg.Output; 

oo(:,ot3x)=oo(:,ot3x)*ot3xcal; oo(:,ot3y)=oo(:,ot3y)*ot3ycal; 

set(figure, 'Name', 'Newport Box 0T3', 'NumberTitle', 'off'); 

plot(oo(:,ot3x),oo(:,ot3y)),grid,zoom 

axis([-5 5 -5 5]) 

xlabel('distance, mm') 

ylabel('distance, mm') 

legend('0T3') 

title('Newport calibration box for OT 3') 

saveas(gcf,'C:\Watkins\MATLAB Files\output figs\Newport calibration box for OT 3.fig') 


% return constants for newport zero on 0T2 


n_y_zero=n_y_zero_ot2; 
n_x_zero=n_x_zero_ot2; 
set (tg,n_x_zero_par,n_x_zero); 
set(tg,n_y_zero_par,n_y_zero); 

newxl=newx2c; newyl=newy2c; n_test_x=const; n_test_y=const; 


save gains calculated 


save mgains b_x_zero b_y_zero bak_xl bak_yl bak_x2 bak_y2 bak_x3 bak_y3 otlxcal 
otlycal... 

ot2xcal ot2ycal ot3xcal ot3ycal newx3c newy3c newx2c newy2c n_x_zero_ot2 
n_x_zero_ot2... 

n_y_zero_ot2 n_y_zero_ot3 finval NGmy NGmx BGmy BGmx 

close_system('calibration_control'); % close the Simulink System 

-tg 


% Clear variables if desired 


reply2 = input('Enter y to clear all variables. Press Enter to continue without 
clearing ','s'); 
if ~isempty(reply2) 
clear 

end 


% Display Results 


disp(['Baker offset, theta X for 0T1 
disp(['Baker offset, theta Y for 0T1 
disp(['Baker theta X gain for 0T1 


',num2str(b_x_zero), ' V']) 

',num2str(b_y_zero), ' V']) 

',num2str(bak_xl), ' V/mm']) 
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disp(['Baker theta Y gain for 0T1 ',num2str(bak_yl),' V/mm']) 

disp(['Baker theta X gain for 0T2 ',num2str(bak_x2),' V/mm']) 

disp(['Baker theta Y gain for 0T2 ',num2str(bak_y2),' V/mm']) 

disp(['Baker theta X gain for 0T3 ',num2str(bak_x3),' V/mm']) 

disp(['Baker theta Y gain for 0T3 ',num2str(bak_y3),' V/mm']) 

disp(['Baker cal const theta Y ',num2str(BGmy),' mRad/V']) 

disp(['Baker cal const theta X ',num2str(BGmx),' mRad/V']) 

disp(['OTl calibration const,X axis ',num2str(otlxcal),' mm/V']) 

disp(['OTl calibration const,Y axis ',num2str(otlycal),' mm/V']) 

disp(['0T2 calibration const,X axis ',num2str(ot2xcal),' mm/V']) 

disp(['0T2 calibration const,Y axis ',num2str(ot2ycal),' mm/V']) 

disp(['0T3 calibration const,X axis ',num2str(ot3xcal),' mm/V']) 

disp(['0T3 calibration const,Y axis ',num2str(ot3ycal),' mm/V']) 

disp(['Newport theta X gain for 0T3, ',num2str(newx3c),' V/mm']) 

disp(['Newport theta Y gain for 0T3, ',num2str(newy3c),' V/mm']) 

disp(['Newport theta X gain for 0T2, ',num2str(newx2c),' V/mm']) 

disp(['Newport theta Y gain for 0T2, ',num2str(newy2c),' V/mm']) 

disp(['Newport cal const theta Y ',num2str(NGmy),' mRad/V']) 

disp(['Newport cal const theta X ',num2str(NGmx),' mRad/V']) 

cd(curr_dir); 

% Clear variables no longer needed 

clear reply replyl reply2 n_test_x_par n_test_y_par b_x_zero_par b_y_zero_par 
n_x_zero_par n_y_zero_par 

clear x_test_par y_test_par bak_x_par bak_y_par newx_par newy_par avgx avgy nlymin 
nlymax 

clear nlxmin nlxmax n3ymax n3ymin n2ymax n2ymin n3xmax n3xmin n2xmax n2xmin otcal calc 
stepsizeyl finall initl 

clear stepsizey2 final2 init2 stepsizey3 final3 init3 stepsizex2 stepsizex3 stepsizexl 
clear valx valy xcnt ycnt valx_sig valy_sig fct fctl fct2 xdist ydist 
format loose 


EOF 
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SAMPLE OUTPUT OF NEWPORT MIRROR CALIBRATION 


I^pcrt 0.5V stqo respcrse,tteta X OT2 axl OT3 


— OT2X 

— OT2Y 

— OT3X 
OT3Y 



ll I_I_I_l_t_t_I 

0 02 0.4 06 Q8 1 1.2 1.4 


1.6 


time 


Q5V step response,theta Y, OT2ardOT3 



Figure 117 X and Y axis Calibration Step Response 

hbAp3trarrpresponse,thetaX OI2andOT3 Nswpot rarrprespcnse.thetaY, OT2ardOT3 




Figure 118 X and Y axis Calibration Ramp Response 
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APPENDIX D: C-CODE FOR ADAPTIVE LATTICE 

ALGORITHMS 


These C-coded algorithms were used in the S-function builder block from 
Simulink. 

A. GRADIENT ADAPTIVE LATTICE 

/* GAL algorithm from Haykin, Table 12.1 */ 

/*nns is the number of stages of the Lattice 
beta is a small constant, less than 0.01 

mu is the adaptation step size, a small number less than 1 
inputs are uO outputs are yO*/ 

int i,n,m; 

{ 

n = *nns; 
m = n + 1; 


/*b0*/ 
/*f0*/ 


yO[0] = uO[7*m]; 
yO[5*m] = uO[7*m]; 


/*y0*/ 

/*ub0*/ 


y0[4*m] = uO[3*m]*u0[7*m]; 
y0[6*m] = *del + uO[7*m]*u0[7*m]; 


/*hl0*/ 


yO[3*m] = uO[3*m] + (*mu)/yO[6*m]*u0[7*m]*u0[7*m+l] ; 


for(i =1; i < m; ++i){ 


/*Em*/ 

/*fm*/ 

/*bm*/ 

/*km*/ 

/*ym*/ 

/*ubm*/ 

/*hml*/ 


yO[2*m+i-l] = (*beta)*u0[2*m+i-l]+(1.0-(*beta))* 

(yO[5*m+i-l]*y0[5*m+i-l]+y0[i—1]*y0[i-1]) ; 

y0[5*m+i] = y0[5*m+i-l] + uO[m+i]*u0[i-1]; 

y0[i] = u0[i-l] + uO[m+i]*y0[5*m+i-l]; 

y0[m+i] = uO[m+i]+(*mu)/yO[2*m+i-l]*(yO[5*m+i-l]* 
yO[i]+u0[i-1]*y0[5*m+i]); 

y0[4*m+i] = y0[4*m+i-l] + yO[3*m+i]*y0[i]; 
y0[6*m+i] = y0[6*m+i-l] + yO[i]*y0[i]; 

y0[3*m+i] = u0[3*m+i] + (*mu)/yO[6*m+i]*y0[i]*u0[7*m+i] ; 

} } 
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B. 


GRADIENT BASED LATTICE 


/* Gradient Based Lattice derived from [Chen, 2001] 
inputs uO, outputs yO 

nst is the number of Lattice stages, lambda is the step size 
int i,n,j; 

double delel[40] , delr[40],delrl[40],ke[40],kel[40],kr[40], 
krl[40],e[40],el[40]; 

double r[40],rl[40],r2[40],ket[40],ketl[40],epsn[40],epsnl[40] ; 
double rho[40],rhol[40],et[40],etl[40],dele[40],yref1; 


n = *nst; 

/*- 

develop input vector 


for (j 


0; j < (n+1); ++j){ 

dele [ j] = uO[j]; 
ke[j] = uO[j+2*(n+1)]; 
e[j] = uO[ j + 4* (n+1)]; 
rl[j] = uO[j + 6*(n+1)]; 
et[j] = uO[ j + 8*(n+1)]; 
rho[j] = uO[j+10*(n+1)]; 


delr[j] = u0[j+n+l]; 
kr[j] = uO[j + 3* (n+1) ] ; 
r [ j] = uO[j + 5* (n+1)] ; 
ket[j] = uO[j + 7* (n+1)]; 
epsn[j] = uO[j+9*(n+1)]; 
} 


■k 


for (j = 0; j < (n+1); ++j) { 

epsnl[j] = epsn[j]; rhol[j] = rho[j]; 

delel[j] = dele[j]; delrl[j] = delr[j]; 

kel[j] = ke[j]; krl[j] = kr[j]; ketl[j] = ket[j]; 

r2[j] = rl [ j]; rl [ j] = r [ j ] ; 

el[j] = e [ j]; etl[j] = et [ j]; } 

r[0] = yhat[0]; e[0] = r[0]; 

et[0] = xhat[0]; 

yref1 = uO[11* (n+1)]; 

/*- 

Algorithm Calculations 

- */ 

for (i=l; i<(n+l); ++i) { 


dele[i]=*lambda*delel[i]+rl[i — 1]*rl[i — 1 ] ; 
delr[i]=*lambda*delrl[i]+e[i-1]*e[i-1]; 

ke[i]=kel[i]+1/delel[i]*r2[i-1]*el [i]; 
kr[i]=krl[i]+1/delrl[i]*el[i-1]*rl[i]; 


e[i] =e[i-1]-rl[i-1]*ke [i]; 
r[i]=rl[i-1]-e[i-1]*kr[i]; 

ket[i]=ketl[i]+l/dele[i]*rl[i-1]*etl[i]; 
et[i]= et[i-1]-r[i-1]*ket [i]; } 


epsn[n] = 0; rho[n] 


0 ; 
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for (j=n-l; j>0 ; —j) { 

epsn[j] = epsn[j+1]-(rho[j+1]+yref[0]*ket[j+1])*kr[j]; 
rho[j] = rhol[j+1]-epsnl[j+1]*ke[j]+yrefl*ket[j+1];} 

ux[0] = -1* (epsn[1]+rho[1]+yref[0]*ket[1]); 


/*- 

/* develop output vector 


for(j = 0; j < (n+1); ++j){ 

y0[j] = dele [ j]; y0[j+n+l] = delr[j]; 

yO[j+2*(n+1)] = ke [ j]; yO[j+3*(n+1)] = 

yO[j + 4*(n+1)] = e [ j]; yO[j + 5*(n+1)] = 

yO[j + 6*(n+1)] = r1 [ j]; yO[j + 7*(n+1)] = 

yO[j + 8*(n+1)] = et [ j]; yO[j+9*(n+1)] = 

yO[j + 10*(n+1)] = rho [ j]; } 

yO[11*(n+1)] = yref [0]; 


*/ 


/* 


kr [ j ] ; 
r [ j] ; 
ket [ j ] ; 
epsn[j] ; 
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APPENDIX E: CALCULATION OF KALMAN AND LQR 

PARAMETERS 


% this program constructs SS rep of 2 axis, 2 input 
% control mirror system 
% States: Vpy, Theta x. Theta-dot x 

% Vpx, Theta y, theta-dot y 

% Gm,Gmd - Mirror Sensor gain, rad/V 
% z,zd - damping; fn,fnd - natural freq, Hz 
% wn,wnd - natural freq, rad/s 

% T,Td - Time delay for mirror/detector response, sec 
"6 

% Control mirror fundamentals 
"6 

T2 = 1/I5e3; % time constant for PSD/input/output system 

o,_ _ __ _ _ _ _ _ 

o 

'natural frequency' of mirrors and mirror damping 


fn2x = 900; 

f n2y 

= 827; z2 = 

0.9; 




wn2x=fn2x*2* 

pi; 

wn2y=fn2y*2 

*pi; 

% convert to rads 

Gm=(52.4e-3) 

720; 

"6 

mirror voltage 

gain. 

rad/V 


Gp=20/10e-3; 


"6 

Position 

Sensor Gain, 

V/m. 


Dm=1245; 


o, 

o 

Distance 

from 

control 

mirror 

to 1 

Dd=0.681; 


"6 

Distance 

from 

disturbance mirror 

Dt=l.5; 


"6 

Distance 

from 

control 

to OT2, 

m . 

Ts=0.0005; 


o, 

o 

Sample time, sec 




Cross-co 

upling terms 




Axl = 

-le-2*wn2x A 2; 




Ax2 = 

-Ie-2*2*wn2x*z2; 




Ayl = 

-le-2*wn2y A 2; 




Ay2 = 

-Ie-2*2*wn2y*z2; 



SS 

Matri 

ces for FSM 



-[ o 


1 

0 

0; 

-wn 

2x A 2 

-2*wn2x*z2 

Axl 

Ax2; 

0 


0 

0 

l; 

Ayl 


Ay2 

-wn2y A 2 

-2*wn2y*z2]; 

;=[0 


0; 



Gm* 

wn2x A 2 

0; 



0 


0; 



0 


Gm*wn2y A 2 ] ; 



: zeros 

(2,4) ; 

C (1,1)=2 *Dm; C(2, 

3)=-C(1,1) 

t 


D=zeros(2,4); 
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Noise Matrices for Kalman Filter to estimate FSM velocity 


Qn 

= le-3*ones (2,2); 

Nn = le-3*ones(2,2); Rn = 2e-l*eye(2,2); 

Gn 

= le-2*ones (4,2) ; 

Gn(2,1) =le-l; Gn (4,2)=Gn(2,1) ; 


Calculate Kalman 

Gains 

L = 

: LQE (A, Gn, C, Qn, Rn, 

Nn) 


Calculate Kalman 

system and convert to discrete 


Ak = (A-L*C); 

Bk = [B L] ; 

Ck = eye(4,4) ; 

Dk = zeros (4,4); 

kalsysc = ss(Ak,Bk,Ck,Dk); 

kalsysd = c2d(kalsysc,Ts); 

ak = kalsysd.a; bk = kalsysd.b; 

ck = kalsysd.c; dk = kalsysd.d; 


%values for geometry considerations of mirrors 



gty = 

0.405; 

gtx : 

= 0.535; 

o 

"o 

State 

Space 

Model for 

Mirror-PSD system 


Asys= 

[ -1/T2 
0 
0 
0 
0 
0 


gty*2*Dm/ot2ycal/T2 0 


0 

-wn2x A 2 

0 

0 

Ayl 


0 


0 


0 ; 


-2*wn2x*z2 

0 

0 

Ay2 


0 0 0 ; 

0 Axl Ax2; 

-1/T2 gtx*2*Dm/ot2xcal/T2 0; 
0 0 1 ; 


0 


-wn2y A 2 


-2*wn2y*z2] 


Bsys=[0 
0 


Gm*wn2x A 2 

0 

0 

0 


0; 

0 ; 

0 ; 

0 ; 

0 ; 

Gm*wn2y A 2]; 


Gsys = zeros(6,2); Gsys(l,l)=l; Gsys(4,2)=l; 


Convert output to mm, ot2xcal/ycal are from calibration program 


Csys = eye(6,6); Dsys = 0; Csys(4,4)=-l; 

Csys2 = zeros(2,6); Csys2(l,l) = ot2ycal; Csys2(2,4)=-ot2xcal; 


Modify SS system to include integrator (after Ogata) 


Cint = zeros (2,6); Cint(l,l) = ot2xcal; Cint(2,4) = ot2ycal; 
Aint = [Asys zeros(6,2);-Cint zeros (2,2)] 

Bint = [Bsys;zeros (2,2)] 
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o\° o\° o\° o\° o\° o\° o\° o\° 


Weighting Matrices for LQRD function 


Qa=1000; 

Q=eye(8,8); Q(l,l)=Qa; Q(4,4)=Qa; 
R=0.l*eye(2,2); 

K = LQRD (Aint, Bint, Q, R, Ts) ; 


Convert SS continuous system to discrete 


dasysc2 = ss(Asys,Bsys,Csys2,zeros (2,2)); 

dasysd2 = c2d(dasysc2,Ts); 

am2=dasysd2.a; 

bm2=dasysd2.b; 

cm2=dasysd2.c; 

dm2=dasysd2.d; 


EOF 
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APPENDIX F: X AXIS RESULTS FOR MULTIPLE FREQUENCY 

PLUS NOISE CASES 


Stile Platform 2 freq PSD Rd - X axis 



Stile Rifom 2 freq PSD Rd-X axis QjTtocotrdles 



Figure 119 PSD of X Axis Response for Stationary Platform, 2 Freq. Disturbance 


StileRatfam+ rase 2freqPSDRd-Xads 



Stile Ratfam+rase 2ft=qPSDRd-X3ds Qrrtocatrdleis 



Figure 120 


PSD of X axis Response for Stationary Platform 2 Freq. plus noise 

Disturbance 
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Power Spectral Density, dB/Hz 


Mbairp Rdfam2frEq PSDHd-Xaas 


Vifcrdirg Rdfarm2 freq PSD Rd-X ads Grrtoccrtrdlas 



— rPXLMB'LGR 

— Q^L/LCR 



0 53 87 100 153 2D 293 

fraqjsncy, hk 


Figure 121 PSD of X Axis Response for Vibrating Platform, 2 Freq. Disturbance 

MbairgRdfarn+rrise 2fraqPSDRd-Xads MbdirgRdfam+nrise 2freqPSDRd-Xads Garboccrtrdlas 




Figure 122 


PSD of X axis Response for Vibrating Platform 2 Freq. plus noise 

Disturbance 
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